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Abstract 

Production of events with hadronic and leptonic final states has been measured in 
e^e^ collisions at centre-of-mass energies of 130-172 GeV, using the OPAL detector at 
LEP. Cross-sections and leptonic forward-backward asymmetries are presented, both in- 
cluding and excluding the dominant production of radiative Z7 events, and compared to 
Standard Model expectations. The ratio i?b of the cross-section for bb production to the 
hadronic cross-section has been measured. In a model-independent fit to the Z lineshape, 
the data have been used to obtain an improved precision on the measurement of 7-Z 
interference. The energy dependence of Oem has been investigated. The measurements 
have also been used to obtain limits on extensions of the Standard Model described by 
effective four-fermion contact interactions, to search for t-channel contributions from new 
massive particles and to place limits on chargino pair production with subsequent decay 
of the chargino into a light gluino and a quark pair. 
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1 Introduction 



Fermion-pair production in e"'"e~ collisions is one of the basic processes of the Standard Model, 
and deviations of measured cross-sections from the predicted values could be a first indication 
of new physics beyond the Standard Model. Measurements up to 161 GeV centre-of-mass 
energy |]1],0 have shown no significant deviations from Standard Model expectations. In this 
paper we present new measurements of hadronic and leptonic final states in e+e" collisions 
at a centre-of-mass energy ^/s of 172 GeV, and improved results for the same final states at 
130, 136, and 161 GeV, using the OPAL detector at LEP. Cross-sections have been measured 
for hadronic, bb, e"'"e~, f^~^n~, and r+r" final states, together with the forward-backward 
asymmetries for the leptonic final states. We present values both including and excluding the 
production of radiative Z7 events. In general, we define a 'non-radiative' sample as events with 
s' /s > 0.8, whereas 'inclusive' measurements are corrected to s'/s > 0.01, where y/s' is defined 
as the centre-of-mass energy of the e"'"e~ system after initial-state radiation. 

In these analyses, we have introduced a well-defined treatment of the interference between 
initial- and final-state photon radiation, and an improved method of taking account of the 
contributions from four-fermion production. While both of these effects are small {0{1%)) 
compared with the statistical precision of the current data, they will become significant with 
the increased luminosity expected at LEP in the future, especially when combining results 
with other experiments [Q. We have reanalysed our data at 130-136 GeV and 161 GeV 
using the same treatment of interference and four-fermion effects, in order to provide a uniform 
sample of measurements for comparison with Standard Model predictions. 

The revised results at 130-136 GeV also benefit from several improvements to the analysis. 
In particular, we benefit from an improved understanding of the background in the inclusive 
hadronic samples arising from two-photon events. The separation of 'non-radiative' hadronic 
events has been improved. The main changes to the lepton analyses include increased efficiency 
for the selection of tau pair events, and the use of a Monte Carlo generator with multiple photon 
emission for simulating the e"'"e~ — >■ e+e~(?T,7) process instead of one containing only single 
photon production. There have also been improvements to the detector calibration, which 
particularly benefit the measurement of R\^, the ratio of the cross-section for bb production to 
the hadronic cross-section. Most of these improvements are already included in the published 
results at 161 GeV 0. 

As has been shown previously ||l|, @], the comparable size of the photon exchange and Z 
exchange amplitudes at these centre-of-mass energies allows constraints to be placed on the 
size of the interference terms between them. In this paper we improve our previous constraints 
by including the data at 172 GeV. In an alternative treatment, we assume the Standard Model 
form of the amplitudes and use the data to investigate the energy dependence of the electro- 
magnetic couphng constant, aem- We have also used the data to search for evidence for physics 
beyond the Standard Model. Firstly we do this within a general framework in which possible 
contributions from extensions of the Standard Model are described by an effective four-fermion 
contact interaction. This analysis is essentially the same as those performed in references 
but the inclusion of data at 172 GeV centre-of-mass energy gives significant improvements to 
the limits presented there. In the previous analysis of the hadronic cross-section we assumed 
the contact interaction was flavour-blind; here we extend the study to include the case where 
the new physics couples exclusively to one up-type quark or one down-type quark. In a second, 
more specific analysis, we set limits on the coupling strength of a new heavy particle which 
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might be exchanged in t-channel production of hadronic final states. Such a particle could be 
a squark, the supersymmetric partner of a quark, in theories where i?-parity is violated, or a 
leptoquark, which is predicted in many theories which connect the quark and lepton sector of 
the Standard Model. In this analysis we assume that the new physics involves only one isomul- 
tiplet of heavy particles coupling with defined helicity. These studies are of topical interest in 
view of the indication of an anomaly at large momentum transfers in e"'"p collisions reported 
by the HERA experiments Contact interactions or production of a heavy particle have 
both been suggested as possible explanations Finally, we place limits on gaugino pair 

production with subsequent decay of the chargino or neutralino into a light gluino and a quark 
pair in supersymmetric extensions to the Standard Model. 

The paper is organized as follows. In section ^ we describe Monte Carlo simulations, the 
treatment of interference effects between initial- and final-state radiation and of the contri- 
butions from four-fermion final states. In section ^ detailed descriptions of the luminosity 
measurement and the analysis of hadronic events, of each lepton channel and of the measure- 
ment of -Rb are given. In section | we compare measured cross-section and asymmetry values 
with Standard Model predictions, and use them to place constraints 7-Z interference and aem- 
Finally, in section |^ we use our measurements to place limits on extensions of the Standard 
Model. 



2 Theoretical Considerations and Simulation 



2.1 Monte Carlo Simulations 



The estimation of efficiencies and background processes makes extensive use of Monte Carlo sim- 
ulations of many different final states. For studies of e^e~ hadrons we used the PYTHIAS. 7 
program with input parameters that have been optimized by a study of global event shape 
variables and particle production rates in Z decay data |TD|. For e+e~ e+e~ we used the 
BHWIDE Monte Carlo program, and for e+e" fi^fi' and e+e" t+t' the KORALZ4.0 
program [|T^. Four-fermion events were modelled with the grc4f |T^, FERMISV and EX- 



CALIBUR generators, with PYTHIA used to check the separate contributions from WW 
and Wez/ diagrams. Two-photon background processes with hadronic final states were simu- 
lated using PYTHIA and PHOJET ]16[ at low Q\ 
with the 'perimiss' option 



At high the TWO GEN program 
was found to give the best description of data; PYTHIA and 



HERWIG ||T9[ were also used for comparison. The Vermaseren generator was used to sim- 
ulate purely leptonic final states in two-photon processes. The e"'"e~ — >■ 77 background in the 
e"'"e~ final state was modelled with the RADCOR program, while the contribution from 



e"'"e~7 where the photon and one of the charged particles are inside the detector acceptance was 
modelled with TEEGG All samples were processed through the OPAL detector simulation 
and reconstructed as for real data. 



program 



23 



For the measurement of the luminosity, the 
cross-section for small-angle Bhabha scattering was calculated using the Monte Carlo program 
BHLUMI , using generated events processed through a simulation program for the forward 
calorimetry. 
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2.2 Initial-final State Photon Interference 



A feature of e^e^ collision data at energies well above the Z mass is a tendency for radiative 
return to the Z. If one or more initial-state radiation photons are emitted which reduce the 
effective centre-of-mass energy of the subsequent e"'"e~ collision to the region of the Z 
resonance, the cross-section is greatly enhanced. In order to test the Standard Model at the 
highest possible energies, we separate clearly radiative events from those with i/i' ~ ^/s using 
methods similar to those in previous analyses p| . In this separation, y/s' is defined as the centre- 
of-mass energy of the e"'"e~ system after initial-state radiation. The existence of interference 
between initial- and final-state radiation means that there is an ambiguity in this definition. 
The Monte Carlo generators used to determine experimental efficiencies and acceptances do 
not include interference between initial- and final-state radiation, but these programs are used 
to correct the data, which do include interference. Therefore further corrections have to be 
applied to the data before measurements can be compared with theoretical predictions. 

For Standard Model predictions (for all channels except e^e~ , which is described below) we 



use the ZFITTER program, which has an option either to enable or to disable interference 
between initial- and final-state radiation. We choose to use the option with interference disabled 
for our comparisons, and correct our measurements to account for this as explained below. This 
choice has the advantage of making the definition of s' unambiguous, and is more suitable for 
interpreting the measurements in terms of theoretical parameters. For example, the S-matrix 



ansatz used to fit the data, described in section 4.1, is unsuitable when the non-resonant part 



of the interference between initial- and final-state radiation contributes PB . 



To determine corrections to the measured cross-sections, we define a differential 'interference 
cross-section' (d^aipsR / dm^j dcos6') as the difference between the differential cross-section 
including initial- final state interference, (d^ajnt / dm^j dcos6'), and that excluding interference, 
(d^cTnoint / dm^j dcos6'), as calculated by ZFITTER using the appropriate fiag settings^. The 
differential interference cross-section may be either positive or negative, depending on the values 
of the cosine of the angle 6 between the fermion and the electron beam direction, and the 
invariant mass of the fermion pair rrifj. We then estimate the fraction of this cross-section 
accepted by our selection cuts by assuming that, as a function of cos 6' and m^j, its selection 
efficiency 

eiFSR(cos6', rriff) = enoint(cos6', nifj), (1) 

where enoint has been determined from Monte Carlo events which do not include interference. 
The corrected cross-section acorr is obtained from the measured cross-section after background 
subtraction and efficiency correction cxmeas as: 



An- An- An f eiFsnfcQS 6^, mf-r) ■'^ "^/^^ a drUf-f 

dcos^ dcos^ dcose Jenoint{cose,mj-^)-£^^f^dmf/ 



In practice the integrals were evaluated in appropriate bins of cos 6 and m^-^. As the accepted 
cross-section is estimated as a function of cos 6', the method is easily applied to total cross- 
sections, angular distributions or asymmetry measurements. 



^Cross-sections including interference are obained by setting INTF=1, those excluding interference by setting 
INTF=0. For hadrons, we also set INCL=0 to enable the INTF flag. 



6 



Interference Corrections (s'/s > 0.8) 




130.25 GeV 


136.22 GeV 


161.34 GeV 


172.12 GeV 


AcThad/o-SM (%) 
Acr^j^/(TSM (%) 
Acr^^/asM (%) 


+ 1.0 ±0.3 ±0.4 
-1.6 ±0.5 
-1.3 ±0.4 


±1.2 ±0.4 ±0.5 
-1.8 ±0.6 
-1.4 ±0.4 


±1.3 ±0.4 ±0.6 
-1.7 ±0.5 
-1.4 ±0.4 


±1.2 ±0.3 ±0.5 
-1.8 ±0.5 
-1.4 ±0.4 



Table 1: Corrections Aa, which have to be applied to the measured non-radiative cross-sections 
in order to remove the contribution from initial-final state interference. They are expressed as 
a percentage fraction of the expected Standard Model cross-section. The first error refiects the 
uncertainty from modelling the selection efficiency for the interference cross-section, the second 
one our estimate of possible additional QCD corrections for the hadrons. 



The systematic error on this procedure was assessed by repeating the estimate modifying 
the assumption of eq. ^ to 

, a ^ enomt(cos6',mff) ± enoint(cos6', v^i) 
eiFSR(cos^,mff) = ^-^ , (3) 

i.e. for each bin of cos 9 and m^j the average of the efficiency in that bin and the efficiency in the 
bin including m^j = ^/s for the same cos 6 range was used. This was motivated as follows. The 
efficiencies enoint used for the interference correction are an average over events with initial-state 
radiation and events with final-state radiation, and are a good approximation to the true eipsR 
if these efficiencies are similar. For large m^j this is the case, but for small m^j the efficiency 
for the relatively rare events with final-state radiation may be significantly higher than that for 
events with initial-state radiation. To account for this, the error on eipsR is taken as half the 
difference between the average efficiency enoint (cos 6*, m^j), and the largest possible efficiency, 
enoint(cos 6*, ^/s) at a givcu value of cos6'. 

For the hadrons there is an additional uncertainty due to QCD effects. We have taken this 
additional uncertainty to be 100% of the correction without an s' cut. The basic assumption 
here is that the near-cancellation between virtual (box) and real interference effects without cuts 



(Kinoshita-Lee-Nauenberg theorem |27|) is not completely destroyed by large QCD corrections 
to both. This has been proven for pure final-state radiation but not yet for initial- final 
state interference. In the absence of a theoretical calculation we allow the asymptotic value to 
change by 100%. 

The corrections to inclusive {s'/s > 0.01) cross-sections are small, refiecting the Kinoshita- 
Lee-Nauenberg cancellation. They typically amount to (— 0.36±0.04)% for muon pairs, (— 0.5± 
0.1)% for tau pairs and (±0.1±0.1)% for hadrons, where the statistical errors are small com- 
pared to the systematic errors, derived as described above. For non-radiative events {s'/s > 0.8) 
the corrections are rather larger, and are given in detail in table |l|. The differences between 
the muon and tau corrections reffect the different acceptance cuts in cos 6* used in the event 
selection; the hadron corrections are of opposite sign from those of the leptons because of the 
quark charges. The corrections change the lepton asymmetry values by typically —0.006 ±0.001 
for s'/s > 0.01 and —0.015 ± 0.005 for s'/s > 0.8. All corrections depend only very weakly on 

VS. 

We have checked the results of the above correction procedure by comparing them to an 
independent estimate using the KORALZ [T^ Monte Carlo generator. We generated two sam- 
ples of muon pair events at 171 GeV, and subjected them to the full detector simulation. 



7 



reconstruction and event selection procedures. Both samples were generated with only single 
photon emission. In the first sample there was no interference between initial- and final-state 
radiation, while in the second sample interference was enabled. The differences between the 
observed cross-sections and asymmetries agreed with the estimates from ZFITTER described 
above within one standard deviation, for both the inclusive {s'/s > 0.01) and non-radiative 
{s'/ s > 0.8) cases. 

The above correction procedure has been applied to all cross-sections, asymmetry measure- 
ments and angular distributions, except for those for the e^e~ final state. In this case, we do 
not use a cut on s' so there is no ambiguity in its definition. Both the Standard Model calcu- 
lations and the Monte Carlo program used to calculate efficiency and acceptance corrections 
include interference between initial- and final-state radiation. Results for e"'"e~ are therefore 
presented including such effects. 

2.3 Four-fermion Effects 

Contributions from four-fermion production fff'f to the process e"'"e~ ff pose non-trivial 
problems both experimentally and theoretically. While four-fermion final states arising from 
the 'two-photon' (multiperipheral) diagrams, for example, can be considered background to two- 
fermion production, those arising from the emission of low mass ff pairs in s-channel diagrams 
may in some circumstances be considered signal, in the same way as is emission of photons. A 
clean separation is not possible because of interference between diagrams contributing to the 
same final state. 

The correct experimental treatment of the four-fermion contribution depends on whether 
or not the theoretical calculation with which the experimental measurement is to be compared 
includes emission of fermion pairs. For example, ALIBABA does not include such emission. 
By default ZFITTER includes initial-state pair emission via virtual photons, although this can 
be disabled. However, pair emission via virtual Z bosons is not included. By comparing 
the predictions of ZFITTER with and without pair emission^, we estimate that the effect of 
including it increases the cross-sections for s'/s > 0.01 by about 1% and decreases those for 
s'/s > 0.8 by about 0.1% at the energies considered here. Similar values are obtained for 
hadrons and lepton pairs. Final-state pair radiation is not explicitly treated in ZFITTER. The 
dominant part of its (very small) effect on the cross-section is covered in the inclusive treatment 
of final state radiation^. For corrections to the selection efficiency, however, both initial- and 
final-state pair radiation have to be considered, as described below. 

None of the theoretical calculations to which we compare our data has an option to separate 
real from virtual fermion pair effects^. Therefore two- and four-fermion events have to be treated 
together in the data analysis. Considering all four-fermion events as background, for example, 
would not account for the virtual vertex corrections, which can be even larger than the effects 
of real pair emission. Therefore some four-fermion events always have to be excluded from 
background estimates. 

^We have modified the ZFITTER code so that the s' cut acts on fermion pairs as well as photons, since by 
default hard pair emission leading to s' < 0.5s is not included. 

^Inclusive treatment of final-state radiation is obtained by setting the flag FINR=:0. 

"^ZFITTER with the flag F0T2=2, like ALIBABA, includes neither real pair emission nor vertex corrections 
involving virtual pairs. The default setting FDT2=3 includes both, summing up beforehand the soft part of real 
pair emission and the vertex corrections. 
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In general, we compare our measurements with ZFITTER predictions including pair emis- 
sion. This means that pair emission via virtual photons from both the initial and final state 
must be included in efficiency calculations, and be excluded from background estimates. In 
order to perform the separation, we ignore interference between s- and t-channel diagrams con- 
tributing to the same four-fermion final state, and generate separate Monte Carlo samples for 
the different diagrams for each final state. For a two-fermion final state ff we then include as 
signal those four-fermion events arising from s-channel processes for which m^j > m^p-, m^,jT < 
70 GeV and rrtj^/s > 0.01 {m'^j/s > 0.8 in the non-radiative case). This kinematic classification 
closely models the desired classification of fff'f in terms of intermediate bosons, in that pairs 
arising from virtual photons are generally included as signal whereas those arising from virtual 
Z bosons are not. All events arising from s-channel processes failing the above cuts, together 
with those arising from the t-channel process (Zee) and two-photon processes are regarded as 
background. Four-fermion processes involving WW or single W production are also background 
in all cases. The overall efficiency, e, is calculated as 

6 = (1 - ^)e,j + ^e,,,jr (4) 

where e^j, Cfff/p- are the efficiencies derived from the two-fermion and four-fermion signal Monte 
Carlo events respectively, CTffpp- is the generated four-fermion cross- sect ion, and cxtot is the total 
cross-section from ZFITTER including pair emission. Using this definition of efficiency, effects 
of cuts on soft pair emission in the four-fermion generator are correctly summed with vertex 
corrections involving virtual pairs. The inclusion of the four-fermion part of the signal produces 
negligible changes to the efficiencies for hadronic events and for lepton pairs with s'/s > 0.8. 
The efficiencies for lepton pairs with s'/s > 0.01 are decreased by about 0.5%. 

The discussion in the above paragraph applies to hadronic, muon pair and tau pair final 
states. In the case of electron pairs, the situation is slightly different. In principle the t-channel 
process with a second fermion pair arising from the conversion of a virtual photon emitted from 
an initial- or final-state electron should be included as signal. As this process is not included 
in any program we use for comparison we simply ignore such events: they are not included as 
background as this would underestimate the cross-section. 



3 Data Analysis 

The OPAL detectoiQ, trigger and data acquisition system are fully described elsewhere p]|-|5^ . 
Data from three separate data-taking periods are used in this analysis: 

• Integrated luminosities of 2.7 pb^^ and 2.6 pb~^ recorded at e~'"e~ centre-of-mass energies 
of 130.25 and 136.22 GeV, respectively, in 1995 (LEP1.5). The energy measurements 
have a common systematic uncertainty of 0.05 GeV 



An integrated luminosity of 10.1 pb ^ recorded at an e+e centre-of-mass energy of 
161.34±0.05 GeV || during 1996. 



^OPAL uses a right-handed coordinate system in which the z axis is along the electron beam direction and 
the X axis is horizontal. The polar angle, 9, is measured with respect to the z axis and the azimuthal angle, (p, 
with respect to the x axis. 
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• An integrated luminosity of approximately 9.3 pb^^ at an e"'"e~ centre-of-mass energy of 
172.3 GeV and 1.0 pb~^ at an energy of 170.3 GeV, recorded during 1996. The data from 
these two energies have been analysed together; the luminosity-weighted mean centre-of- 
mass energy has been determined to be 172.12±0.06 GeV |p5|| . 

3.1 Measurement of the Luminosity 

The integrated luminosity was measured using small-angle Bhabha scattering events, e"'"e~ — >■ 
e"'"e~, recorded in the forward calorimetry. The primary detector is a silicon-tungsten lumi- 
nometer |^ which consists of two finely segmented silicon-tungsten calorimeters placed around 
the beam pipe, symmetrically on the left and right sides on the OPAL detector, 2.4 m from 
the interaction point. Each calorimeter covers angles from the beam between 25 and 59 mrad. 
Bhabha scattering events were selected by requiring a high energy cluster in each end of the 
detector, using asymmetric acceptance cuts. The energy in each calorimeter had to be at least 
half the beam energy, and the average energy had to be at least three quarters of the beam 
energy. The two highest energy clusters were required to be back-to-back in 0, | — 0l| — 7r| < 
200 mrad, where (f)R and (pL are the azimuthal angles of the cluster in the right- and left-hand 
calorimeter respectively. They were also required to be coUinear, by placing a cut on the differ- 
ence between the radial positions, Ai? = \Rr — Rl\ < 25 mm, where Rr and Rl are the radial 
coordinates of the clusters on a plane approximately 7 radiation lengths into the calorimeter. 
This cut, corresponding to an acollinearity angle of about 10.4 mrad, effectively defines the 
acceptance for single-photon radiative events, thus reducing the sensitivity of the measurement 
to the detailed energy response of the calorimeter. The distribution of Ai? for the 172 GeV 
data is shown in figure 0(a). 

For the 130-136 GeV data, the inner and outer radial acceptance cuts delimited a region 
between 31 and 52 mrad on one side of the calorimeter, while for the opposite calorimeter a 
wider zone between 27 and 56 mrad was used. Two luminosity measurements were formed with 
the narrower acceptance on one side or the other side. The final measurement was the average 
of the two and has no first order dependence on beam offsets or tilts. Before data-taking started 
at -(/i=161 GeV, tungsten shields designed to protect the tracking detectors from synchrotron 
radiation were installed around the beam pipe. The shields, 5 mm thick and 334 mm long, 
present roughly 50 radiation lengths to particles originating from the interaction region, almost 
completely absorbing electromagnetically showering particles between 26 and 33 mrad from the 
beam axis. The fiducial regions for accepting Bhabha events for the 161 and 172 GeV data 
were therefore reduced, to between 38 and 52 mrad on one side and between 34 and 56 mrad 
on the opposite side. The distributions of the radial coordinates of the clusters for the 172 GeV 
data are shown in figure ||(b,c). 

The error on the luminosity measurement is dominated by data statistics. For the 130 
and 136 GeV data, the acceptance of the luminometer was reduced at the trigger level by a 
prescaling factor of 16 in order to increase the experimental live time as far as possible, giving 
a statistical error of 0.9% on the combined 130 and 136 GeV data. For the two higher energies 
this prescaling factor was reduced to 2 or 4, and the statistical error amounts to 0.42% (0.43%) 
at 161 (172) GeV. The largest systematic uncertainty arises from theoretical knowledge of 
the cross-section (0.25%), with detector effects amounting to a further 0.20% (0.23%) at 161 
(172) GeV. 
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A second luminosity measurement was provided by the forward detector, a lead-scintillator 
sampling calorimeter covering angles from the beam between 40 and 150 mrad. The selection 



of Bhabha events within the calorimeter acceptance is unchanged from reference p6[ , but the 



acceptance was reduced to the region between 65 and 105 mrad from the beam because of the 
addition of the silicon-tungsten luminometer on the inside front edge of the device. The overall 
acceptance of the calorimeter was measured by normalizing to the precisely known cross-section 
for hadronic events at the Z peak, and applying small corrections derived from Monte Carlo 
simulations to reflect changes in acceptance with centre-of-mass energy. To allow for changes 
in acceptance between years, this normalization was performed separately for 1995 and 1996 
using data recorded at the Z in each year. Knowledge of the hadronic acceptance for the Z data 
is the main source of systematic error in the forward detector luminosity measurement, which 
amounts to 0.8% (1.0%) for the data taken in 1995 (1996). 

The luminosity measured by the forward detector agreed with that measured by the silicon- 
tungsten luminometer to within one standard deviation of the combined error for all data 
samples. For the 130 and 136 GeV data, where the precision of the two measurements was 
similar, the average luminosity was used; the overall error on this average measurement is 
0.7%. At 161 and 172 GeV the silicon-tungsten luminosity was preferred as the more precise; 
the overall error on this measurement amounts to 0.53% (0.55%) at 161 (172) GeV. The errors 
on luminosity are included in the systematic errors on all cross-section measurements presented 
in this paper. Correlations between cross-section measurements arising from errors in the 
luminosity have been taken into account in the interpretation of the results. 



3.2 Hadronic Events 

3.2.1 Inclusive Events {s'/s > 0.01) 

The criteria used to select an inclusive sample of hadronic events with s'/ s > 0.01 were based on 
energy clusters in the electromagnetic calorimeter and the charged track multiplicity. Clusters 
in the barrel region were required to have an energy of at least 100 MeV, and clusters in the 
endcap detectors were required to contain at least two adjacent lead glass blocks and have an 
energy of at least 200 MeV. Tracks were required to have at least 20 measured space points. 
The point of closest approach to the nominal beam axis was found, and required to lie less than 
2 cm in the r-0 plane and less than 40 cm along the beam axis from the nominal interaction 
point. Tracks were also required to have a minimum momentum component transverse to the 
beam direction of 50 MeV. 

The following requirements were used to select hadronic candidates. 

• To reject leptonic final states, events were required to have high multiplicity: at least 
7 electromagnetic clusters and at least 5 tracks. 

• Background from two-photon events was reduced by requiring a total energy deposited 
in the electromagnetic calorimeter of at least 14% of the centre-of-mass energy: R^^s = 
Si^cius/v^ > 0.14, where Edus is the energy of each cluster. 

• Any remaining background from beam-gas and beam-wall interactions was removed, and 
two-photon events further reduced, by requiring an energy balance along the beam direc- 
tion which satisfied -Rbai =| S(_Ecius ■ cos6') | /T^Edns < 0.75, where 9 is the polar angle of 
the cluster. 
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These criteria are identical to those used previously at 161 GeV 0, but the cut on i?bai is 
somewhat looser than that used previously at 130-136 GeV, resulting in a slightly higher 
efficiency for radiative return events. Distributions of i?vis and i?bai for each centre-of-mass 
energy are shown in figure |[ The efficiency of the selection cuts was determined from Monte 
Carlo simulations, and the value for each centre-of-mass energy is given in table 0. From 
comparisons of the data distributions of i?bai and i?vis with Monte Carlo, at these energies 
and at energies around the Z peak (LEPl), we estimate the systematic error on the selection 
efficiency to be 1%. 

Above the W-pair threshold, the largest single contribution to the background arises from 
WW events. No cuts have been applied to reject W-pair events; the expected contribution 
from these to the visible cross-section has been subtracted, and amounts to (2.4±0.2)% at 
161 GeV and (9.6±0.2)% at 172 GeV, where the error arises mainly from the uncertainty in 
the W mass Backgrounds to the inclusive hadron samples at all energies arise from other 
four-fermion events which are not considered part of the signal, in particular two-photon events 
and the channels Zee and Wez/, and tau pairs. These amount to 1.9% at 130 GeV, rising to 
4.1% at 172 GeV. The main uncertainty on this background arises from the two-photon events; 
we assign a 50% error to this contribution, which covers the predictions from all the generators 
discussed in section ETT. 



The numbers of selected events and the resulting cross-sections are shown in table 



3.2.2 Non-radiative events {s'/s > 0.8) 

The effective centre-of-mass energy \/s' of the e"'"e~ collision for hadronic events selected as 
above was estimated as follows. The method is the same as that used in reference 0. Isolated 
photons in the electromagnetic calorimeter were identified, and the remaining tracks, electro- 
magnetic and hadron calorimeter clusters formed into jets using the Durham (fc-r) scheme 
with a jet resolution parameter y^ut = 0.02. If more than four jets were found the number was 
forced to be four. The jet energies and angles were corrected for double counting using the 



algorithm described in reference [^. The jets and observed photons were then subjected to a 
series of kinematic ffis with the constraints of energy and momentum conservation, in which 
zero, one, or two additional photons emitted close to the beam direction were allowed. The fit 
with the lowest number of extra photons which gave an acceptable "w^as chosen. The value 
of \/s' was then computed from the fitted four-momenta of the jets, i.e. excluding photons 
identified in the detector or those close to the beam direction resulting from the ffi. If none 
of the kinematic ffis gave an acceptable was estimated directly from the angles of the 

jets as in reference |jl|. Figure ^ shows distributions at 172 GeV for events with different 
numbers of photons. Note that this algorithm results in s' equal to s for events which give a 
good kinematic ffi with no photon either in the detector or along the beam direction. 

Non-radiative events were selected by demanding s'/ s > 0.8. The numbers of events selected 
at each energy are shown in table 0, together with the corresponding efficiencies and the frac- 
tions of the s'/s > 0.8 sample arising from feedthrough of events with lower s'/s, determined 
from Monte Carlo simulations. 

The estimation of background in the non-radiative samples is less problematic than in the 
inclusive case, because the contribution from two-photon events is tiny. The largest contri- 
bution arises from W-pair events (above the W-pair threshold), and as in the inclusive case 
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Hadrons (s'/s > 0.01) 




130.25 GeV 


136.22 GeV 


161.34 GeV 


172.12 GeV 


Events 


832 


673 


1472 


1368 


Efficiency (%) 


95.8±1.0 


95.2±1.0 


92.3±0.9 


91.2±0.9 


Background (pb) 


5.9±2.1 


5.9±2.1 


8.5±1.7 


15.6±1.5 


CTmeas (pb) 


317±11±5 


264±10±4 


150±4±2 


127±4±2 


C^corr (pb) 


317±11±5 


264±10±4 


150±4±2 


127±4±2 


o"SM (pb) 


330 


273 


150 


125 


Hadrons {s'/s > 0.8) 




130.25 GeV 


136.22 GeV 


161.34 GeV 


172.12 GeV 


Events 


174 


166 


370 


339 


Efficiency (%) 


91.0±0.7 


91.0±0.7 


91.8±0.5 


91.8±0.5 


Feedthrough (%) 


8.9±1.9 


7.9±1.9 


5.2±1.9 


4.8±1.9 


Background (pb) 


1.3±0.1 


1.3±0.1 


2.73±0.22 


6.82±0.26 


0"meas (pb) 


63.5±4.9±1.5 


63.1±5.0±1.5 


35.1±2.0±0.8 


26.7±1.8±0.6 


O"corr (pb) 


64.3±4.9±1.5 


63.8±5.0±1.5 


35.5±2.0±0.8 


27.0±1.8±0.6 


o"SM (pb) 


77.6 


62.9 


33.7 


27.6 



Table 2: Numbers of selected events, efficiencies, backgrounds, feedthrough of events from lower 
s' to the s'/s > 0.8 samples and measured cross-sections for hadronic events. The errors on 
efficiencies, feedthrough and background include Monte Carlo statistics and systematic effects, 
with the latter dominant. The cross-sections labelled cr^eas are the measured values without 
the correction for interference between initial- and final-state radiation, those labelled (Tcorr are 
with this correction. For the inclusive measurements, the results are the same to the quoted 
precision. The ffist error on each measured cross-section is statistical, the second systematic. 
The Standard Model predictions, asm, are from the ZFITTER p5[ program. 



the expected contribution has been subtracted. This amounts to (6.0±0.5)% at 161 GeV and 
(19.8±0.2)% at 172 GeV, where again the dominant error reflects the uncertainty in the W 
mass. Additional small backgrounds arise from four-fermion production and tau pair events. 
The total background at each energy is shown in table ^, together with the final non-radiative 
hadronic cross-sections. 

The main systematic error arises from the modelling of the separation of non-radiative from 
clearly radiative events, and was estimated by comparing eight different methods of separation. 
For example, the algorithm was changed to allow for only a single radiated photon, the photon 
identification algorithm was modified, the hadron calorimeter was removed from the analysis or 
the jet resolution parameter was altered. In each case, the modified algorithm was applied to 
data and Monte Carlo, and the corrected cross-section computed. The changes observed were 
in all cases compatible with statistical fluctuations, but to be conservative the largest change 
(averaged over all beam energies) was taken to define the systematic error, amounting to 2.0%. 
This error is expected to decrease in future with improved data statistics. The error arising from 
the subtraction of W-pair background was investigated by performing an alternative analysis in 
which events identified as W-pairs according to the criteria in reference were rejected. The 
resulting cross-sections after correcting to no interference, 35.4±1.9±0.8 (26.5±1.7±0.6) pb at 
161 (172) GeV are in excellent agreement with those obtained by subtracting the expected 
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W-pair contribution. 

To measure the angular distribution of the primary quark in the hadronic events, we have 
used as an estimator the thrust axis for each event determined from the observed tracks and 
clusters. The angular distribution of the thrust axis was then corrected to the primary quark 
level using bin-by-bin corrections determined from Monte Carlo events. No attempt was made 
to identify the charge in these events, and thus we measured the folded angular distribution. 
The measured values for the s' / s > 0.8 sample are given in table |^. 



Hadrons (s'/s > 0.8) 


cos^^ 


dcr/d cos^^ (pb) 




130.25 GeV 


136.22 GeV 


161.34 GeV 


172.12 GeV 


[0.0,0.1] 


44±13 


35±11 


28.0±5.5 


22.3±5.2 


[0.1,0.2] 


47±12 


57±14 


26.9±5.4 


21.7±5.3 


[0.2,0.3] 


58±15 


44±13 


36.8±6.3 


18.0±4.9 


[0.3,0.4] 


51±13 


58±15 


28.4±5.5 


27.5±5.9 


[0.4,0.5] 


52±14 


38±12 


24.3±5.1 


26.5±5.8 


[0.5,0.6] 


53±14 


53±14 


29.9±5.6 


13.8±4.4 


[0.6,0.7] 


70±16 


114±20 


42.5±6.6 


36.1±6.4 


[0.7,0.8] 


66±15 


57±14 


45.5±6.8 


32.5±5.9 


[0.8,0.9] 


84±17 


69±16 


38.5±6.2 


34.6±6.1 


[0.9,1.0] 


141±31 


122±29 


57.7±10.3 


35.6±8.2 



Table 3: Differential cross-sections for qq production. The values are corrected to no interfer- 
ence between initial- and final-state radiation as described in the text. Errors include statistical 
and systematic effects combined, with the former dominant. 



3.3 Electron Pairs 

The production of electron pairs is dominated by t-channel photon exchange, for which a 
definition of s' as for the other channels is less natural. In addition, the increased probability 
for final-state radiation relative to initial-state radiation renders the separation between initial- 
and final-state photons more difficult. Events with little radiation were therefore selected by a 
cut on 6'acoi, the acollinearity angle between electron and positron. A cut of ^^acoi < 10° roughly 
corresponds to a cut on the effective centre-of-mass energy of s'/s > 0.8, for the s-channel 
contribution. We measure cross-sections for three different acceptance regions, defined in terms 
of the angle of the electron, 9^- , or positron, 9^+ , with respect to the incoming electron direction, 
and the acollinearity angle: 

• A: I cos 6^0- 1 < 0.9, | cos 6*0+ 1 < 0.9, ^^acoi < 170°; this is a loose 'inclusive' measurement; 

• B: I cos6'e-| < 0.7, ^acoi < 10°; this acceptance region is expected to be enriched in the s- 
channel contribution, and is used for asymmetry measurements; in addition, we measure 
the angular distribution for | cos^^c-| < 0.9 and 6'acoi < 10°; 
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• C: I cos6'e-| < 0.96, | cos 6*0+1 < 0.96, 6'acoi < 10°; this 'large acceptance' selection acts as 
a check on the luminosity measurements. 



The selection of electron pair events is identical to previous analyses 0. Events selected 
as electron pairs are required to have at least two and not more than eight clusters in the 
electromagnetic calorimeter, and not more than eight tracks in the central tracking chambers. 
At least two clusters must have an energy exceeding 20% of the beam energy, and the total 
energy deposited in the electromagnetic calorimeter must be at least 50% of the centre-of-mass 
energy. For selections A and B, at least two of the three highest energy clusters must each 
have an associated central detector track. If a cluster has more than one associated track, the 
highest momentum one is chosen. If all three clusters have an associated track, the two highest 
energy clusters are chosen to be the electron and positron. For the large acceptance selection, 
C, no requirement is placed on the association of tracks to clusters, but the requirement on the 
total electromagnetic energy is increased to 70% of the centre-of-mass energy. 

These cuts have a very high efficiency for e"^e~ events while providing excellent rejection of 
backgrounds, which either have high multiplicity or lower energy deposited in the electromag- 
netic calorimeter. The efficiency of the selection cuts, and small acceptance corrections, have 
been determined using Monte Carlo events generated with the BHWIDE []TT| program. These 



are found to be independent of energy over the range considered here. Remaining backgrounds 
arise from t~^t~ events and, in the case of the loose acollinearity cut, also from electron pairs in 
two-photon events and from radiative Bhabha scattering events in which one electron is outside 
the detector acceptance but the photon is within the acceptance. In the case of the large ac- 
ceptance selection, C, which does not require tracks, the main background arises from 77 final 
states. The efficiencies and backgrounds at the three energies are summarized in table ^ In 
figure ^(a,b) we show distributions of total electromagnetic calorimeter energy, after all other 
cuts, for acceptance regions B and C at 172 GeV, showing reasonable agreement between data 
and Monte Carlo. The degraded energy resolution in acceptance region C arises from the in- 
creased amount of material in front of the electromagnetic calorimeter at large | cos^|, where 
the events are concentrated. The acollinearity angle distribution for the inclusive selection. A, 
is shown in figure §(c), and we see good agreement between data and Monte Carlo expectation, 
including the peak corresponding to radiative s-channel return to the Z. 

The numbers of selected events and resulting cross-sections are shown in table ^. The 
following sources of systematic error in the cross-section measurements have been considered. 



• Deficiencies in the simulation of the selection cuts. As shown in figure |5(b), the total 
calorimeter energy distribution is slightly broader in data than Monte Carlo for the large 
acceptance selection, C. The effect of this on the efficiency of this selection has been 
estimated by varying the cut in the range 40% to 75% of the centre-of-mass energy. In 
the other two selections, a more important effect is the efficiency for finding two tracks, 
which has been investigated using events in which only one cluster has an associated 
track. 

• Knowledge of the acceptance correction and how well the edge of the acceptance is mod- 
elled. Because of the steeply falling distribution, any bias in the measurement of 6 has 
a significant effect on the cross-sections, particularly for the large acceptance selection. 
This has been investigated by comparing measurements of 6 made using central detector 
tracks, calorimeter clusters and the outer muon chambers. In addition, in each case the 
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full size of the acceptance correction derived from Monte Carlo has been included as a 
systematic error. 

• Uncertainties in the background contributions. For selections A and B these have been 
assessed by considering the numbers of events failing the total energy cut. Data and Monte 
Carlo are consistent within the statistical precision of 30%. The background in the large 
acceptance selection is almost all from 77 final states, which is much less uncertain. 



The total systematic error in selections A and B amounts to 1.4% and 0.8% respectively, of 
which the largest contribution arises from uncertainty in the track matching efficiency (0.8% 
and 0.5% respectively). In the large acceptance selection, the largest component in the total 
systematic error of 1.1% arises from uncertainty in modelling the edge of the acceptance (0.9%). 



e+e-(A: | cos0| < 0.9, 0acoi < 170°) 




130.25 GeV 


136.22 GeV 


161.34 GeV 


172.12 GeV 


Events 


591 


514 


1587 


1397 


Efficiency (%) 




98.2±1.3 




Background (pb) 


3.7±1.1 


3.4±1.0 


2.3±0.7 


1.9±0.6 


O-mcas (pb) 


220±9±3 


197±9±3 


158±4±2 


135±4±2 


o"SM (pb) 


237 


217 


154 


135 


e+e-(B: | cos6le-| < 0.7, 6>acoi < 10°) 




130.25 GeV 


136.22 GeV 


161.34 GeV 


172.12 GeV 


Events 


112 


98 


285 


246 


Efficiency (%) 




99.2±0.7 




Background (pb) 


0.6±0.2 


0.5±0.2 


0.4±0.1 


0.3±0.1 


CTmcas (pb) 


41.3±4.0±0.5 


37.3±3.8±0.4 


28.1±1.7±0.3 


23.5±1.5±0.2 


c^SM (pb) 


43.1 


39.5 


28.1 


24.7 


e+e-(C: | cos0| < 0.96, 0acoi < 10°) 




130.25 GeV 


136.22 GeV 


161.34 GeV 


172.12 GeV 


Events 


1686 


1542 


4446 


3870 


Efficiency (%) 




98.5±1.1 




Background (pb) 


21.1±2.1 


19.3±1.9 


13.9±1.4 


12.2±1.2 


CTmeas (pb) 


615±16±8 


580±15±8 


434±7±5 


365±6±5 


crsM (pb) 


645 


592 


425 


375 



Table 4: Numbers of selected events, efficiencies, backgrounds and measured cross-sections for 
e"'"e~ events. The efficiencies are effective values combining the efficiency of selection cuts for 
events within the acceptance region and the effect of acceptance corrections. The errors on 
the efficiencies and backgrounds include Monte Carlo statistics and all systematic effects, the 
latter being dominant. The first error on each measured cross-section is statistical, the second 
systematic. The Standard Model predictions, (Tsm, are from the ALIBABA program. 



Unlike all other channels, values for e"'"e include the effect of interference between initial- and 
final-state radiation. 



The measurement of the angular distribution and asymmetry uses the same event selection 
as above, with the further requirement that the two tracks have opposite charge. This extra 
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requirement reduces the efficiency by about 2.5% in the region | cos6\ < 0.9. In addition, to 
reduce the effect of charge misassignment, events with cos^e- < —0.8 must satisfy two extra 
criteria: both electron and positron tracks must have momenta of at least 25% of the beam 
momentum, and there must be only one good track associated with each cluster. These extra 
criteria reduce the overall correction factor to the angular distribution for cos6'e- < —0.8 from 
about 25% to 5%. 

The observed angular distribution of the electron, for events with 6'acoi < 10°, is shown in 
figure §(d). As the variation of the angular distribution with energy is small over the range 
considered here, we have summed data from all energies for this comparison with Monte Carlo 
expectation. The corrected distributions in cos 6 at each energy are given in table ^. Systematic 
errors, arising mainly from uncertainty in the efficiency for finding two tracks with opposite 
charge, amount to 1.2% and are included in the errors in table ^ The forward-backward asym- 
metries for the ^acoi < 10° sample at each energy within the angular range | cos0c-| < 0.7 were 
evaluated by counting the numbers of events in the forward and backward cos 6^- hemispheres. 
The measured values are shown in table |[ Again, the errors are predominantly statistical, 
with small systematic effects arising from charge misassignment, acceptance definition and 
background included in the values given. 

In figure ^(b) we show the distribution of y/s' for the inclusive electron pair events at 
172 GeV. The value of s' for each event was estimated from the polar angles, 6i and 62, of the 
two electrons, assuming massless three-body kinematics to calculate the energy of a possible 
undetected initial-state photon along the beam direction as ■ | sin(6'i -|- ^2) |/(| sin(6'i + 6^2) | + 
sin^i -t- sin ^2)- A similar technique was used to calculate s' for muon pairs and tau pairs. In 
contrast to the other final states, the radiative return peak forms only a very small contribution 
to this channel. 
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Table 5: Differential cross-sections for lepton pair production. The values for e+e^ are for 
^acoi < 10°; those for and r+r" are for s'/s > 0.8 and are corrected to no interference 

between initial- and final-state radiation as described in the text. Errors include statistical and 
systematic effects combined, with the former dominant. 



18 



e+e-(| COS0e-| < 0.7, ^acol < 10°) 





130.25 GeV 


136.22 GeV 


161.34 GeV 


172.12 GeV 




98 


84 


257 


222 


Nb 


12 


13 


17 


17 


A meas 
^FB 


0.79±0.06 


0.73±0.07 


0.88±0.03 


0.86±0.04 


/iSM 
^FB 


0.80 


0.80 


0.81 


0.81 



Table 6: The numbers of forward (Np) and backward (iVs) events and measured asymmetry 
values for electron pairs. The measured asymmetry values include corrections for background 
and efficiency. The errors shown are the combined statistical and systematic errors; in each 
case the statistical error is dominant. The Standard Model predictions, ^|^, are from the 



ALIBABA |29| program. Unlike all other channels, values for e"'"e include the effect of inter- 
ference between initial- and final-state radiation. 



3.4 Muon Pairs 

The selection of muon pair events is essentially identical to previous analyses except that 
the cut on visible energy has been made dependent on the centre-of-mass energy, to reduce loss 
of radiative return events at higher energies. Muon pair events were required to have at least 
two tracks with momentum greater than 6 GeV and |cos^| < 0.95, separated in azimuthal 
angle by more than 320 mrad, and identified as muons. These tracks must have at least 20 
hits in the central tracking chambers and the point of closest approach to the nominal beam 
axis must lie less than 1 cm in the r-0 plane and less than 50 cm along the beam axis from 
the nominal interaction point. To be identified as a muon, a track had to satisfy any of the 
following conditions: 

• At least 2 muon chamber hits associated with the track within A0 = (100 + 100/p) mrad, 
with the momentum p in GeV; 

• At least 4 hadron calorimeter strips associated with the track within = (20 + 100/p) mrad, 
with p in GeV. The average number of strips in layers containing hits had to be less than 

2 to discriminate against hadrons. For | cos^| < 0.65, where tracks traverse all 9 layers of 
strips in the barrel calorimeter, a hit in one of the last 3 layers of strips was required; 

• Momentum p > 15 GeV and the electromagnetic energy associated to the track within 
A(f) < 70 mrad less than 3 GeV. 

If more than one pair of tracks satisfied the above conditions, the pair with the largest total 
momentum was chosen. Background from high multiplicity events was rejected by requiring 
that there be no more than one other track in the event with a transverse momentum greater 
than 0.7 GeV. 

Background from cosmic ray events was removed using the time-of-fiight (TOF) counters 
and vertex cuts. In the barrel region, at least one TOF measurement was required within 
10 ns of that expected for a particle coming from the interaction point. In addition, back-to- 
back pairs of TOF counters were used to reject cosmic rays which had traversed the detector. 
Figure ^(a) shows the distribution of time difference. At, between pairs of back-to-back TOF 
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counters for muon pair events, before applying this cut, clearly showing one peak at the origin 
from muon pairs and a second peak at about 15 ns from cosmic rays. In the forward region, 
for which TOF information was not available, the matching of the central detector tracks to 
the interaction vertex was used in order to remove cosmic ray background. The cosmic ray 
contamination after all cuts is low. There are no events remaining close to the cosmic ray 
rejection cut boundaries in the 130 and 136 GeV samples, and one event remaining in each of 
the 161 and 172 GeV samples. 

Background from two-photon events was rejected by placing a cut on the total visible energy, 
E'vis, defined as the scalar sum of the momenta of the two muons plus the energy of the highest 
energy cluster in the electromagnetic calorimeter: 

Rvis = E,,J^s > 0.5(m|/s) + 0.35. 

The value of this cut is 0.15 below the expected value of i?vis for muon pairs in radiative return 
events where the photon escapes detection, visible as secondary peaks in figure H(b-d). 

The value of s' for each event was estimated from the polar angles of the two muons, as 
described in section |3.3| for electrons. The observed distribution of x/y at 172 GeV is shown 



in figure ^(c). A non-radiative sample of events was selected by requiring s'/s > 0.8. The 
selection efficiencies, and feedthrough of events from lower s'/s into the non-radiative sample, 
were determined from Monte Carlo simulations, and are shown in table |^. 

The residual background in the inclusive sample, of around 4% at 130 GeV increasing to 
11% at 172 GeV, arises mainly from e"'"e~/x"'"/i~ final states, while that in the non-radiative 
sample is predominantly t~^t~ events and amounts to about 5% in total. Total backgrounds 
are shown in table ^ together with the numbers of selected events and resulting cross-section 
measurements. 

Systematic errors on the cross-section measurements, which arise from uncertainties in effi- 
ciency and backgrounds, are small compared to the statistical errors. In all cases, the dominant 
systematic error arises from the uncertainty in the background contamination. 

The observed angular distribution of the fi~ is shown in figure 0(a) for the s'/s > 0.01 
sample and figure |^(b) for the s'/s > 0.8 sample, for all centre-of-mass energies combined. The 
angular distributions at each energy were corrected for efficiency and background, including 
feedthrough of muon pair events from lower s'/s into the non-radiative samples, using Monte 
Carlo events. The corrected angular distributions are shown in table ^. The final values 
have been obtained by averaging the distribution measured using the negative muon with that 
using the positive muon; although this averaging does not reduce the statistical errors on 
the measurements, it is expected to reduce most systematic effects. The forward-backward 
asymmetries at each energy were obtained by counting the numbers of events in the forward 
and backward hemispheres, after correcting for background and efficiency. The data at 130 
and 136 GeV have been combined for the asymmetry measurements. Systematic errors were 
assessed by comparing results obtained using different combinations of tracking and muon 
chambers to measure the muon angles. The total systematic error, including the contribution 
from the correction for interference between initial- and final-state radiation, is below 0.01 in all 
cases, much smaller than the statistical errors. The measured asymmetry values are compared 
with the Standard Model predictions in table & 
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Table 7: Numbers of selected events, efficiencies, backgrounds, feedthrougli of events from lower 
s' to the s'/s > 0.8 samples and measured cross-sections for muon pair events. The errors on 
efficiencies and background include Monte Carlo statistics and systematic effects. The cross- 
sections labelled o"meas are the measured values without the correction for interference between 
initial- and final-state radiation, those labelled (Jcorr are with this correction. In some cases, 
the results are the same to the quoted precision. The ffist error on each measured cross- 
section is statistical, the second systematic. The Standard Model predictions, (Jsm, are from 
the ZFITTER ||25| program. 



3.5 Tau Pairs 



The selection of e+e~ — > r+r~ events is based on that used in previous analyses [|T],3, using 
information from the central tracking detectors and electromagnetic calorimetry to identify 
events with two collimated, low multiplicity jets. However, the cuts have been optimized and 
unified for the different energies, giving an improved efficiency at 130-136 GeV. An inclusive 
sample of events was selected with the following cuts. 



• Hadronic events were rejected by demanding low multiplicity: the number of tracks re- 
constructed in the central tracking detectors had to be at least two and at most six, and 
the sum of the number of tracks and the number of electromagnetic clusters not more 
than 15. 

• The total energy of an event was restricted in order to reject events from e"'"e~ — >■ e"'"e~(7) 
and two-photon processes: the total visible energy in the event, derived from the scalar 
sum of all track momenta plus electromagnetic calorimeter energy, was required to be 
between 0.3y/s and l.ly/s. In addition, the total electromagnetic calorimeter energy was 
required to be less than 0.7y/s and the scalar sum of track momenta less than 0.8-y/s. 
In the endcap region, | cos 6^1 > 0.7, the upper limit on the visible energy was reduced 
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Table 8: The numbers of forward (A^f) and backward (A'^b) events and measured asymme- 
try values for muon and tau pairs. The measured asymmetry values include corrections for 
background and efficiency and are corrected to the full solid angle. The errors shown are the 
combined statistical and systematic errors; in each case the systematic error is less than 0.01. 
The values labelled Ap^^^ are the measured values without the correction for interference be- 
tween initial- and final-state radiation, those labelled A^^^ are with this correction. In some 
cases, the results are the same to the quoted precision. The Standard Model predictions, A^q, 
are from the ZFITTER program. 
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to I.OSa/s because of the less good electron energy resolution. The distribution of total 
visible energy, after all other cuts have been applied, is shown in figure §(a) for all centre- 
of-mass energies combined. 

Background from two-photon events was further reduced by cuts on the missing mo- 
mentum and its direction. The missing momentum in the transverse plane to the beam 
axis, calculated using the electromagnetic calorimeter, was required to exceed 1.5 GeV. 
The polar angle of the missing momentum calculated using tracks only or electromag- 
netic clusters only was required to satisfy | cos6'| < 0.95 and | cos^^| < 0.875 respectively. 
Figure ^(b) shows the distribution of cos^^ of the missing momentum vector calculated 
using electromagnetic clusters after all other cuts have been applied, for all centre-of-mass 
energies combined. 

Vertex and TOF cuts were imposed to remove cosmic ray events, as for events. 
In addition, e"'"e~ fi~^^~ events were removed; these were identified by the criteria 
described in section R]4|, except that the total visible energy was required to exceed 60% 



of the centre-of-mass energy. 

Tau pair events are characterized by a pair of narrow 'jets'. Tracks and electromagnetic 
clusters, each treated as separate particles, were combined in the following way. First the 
highest energy particle in the event was selected and a cone with a half angle of 35° was 
defined around it. The particle with the next highest energy inside the cone was combined 
with the first. The momenta of the combined particles were added and the direction of 
the sum was used to define a new cone, inside which the next highest energy particle was 
again looked for. This procedure was repeated until no more particles were found inside 
the cone. Similarly, starting with the highest energy particle among the remainder, a new 
cone was initiated and treated in the same way. This process continued until finally all 
the particles in the event had been assigned to a cone. 

At least one charged particle was required for each cone, and the sum of the energy in the 
electromagnetic calorimeter and the track momenta in a cone had to be more than 1% of 
the beam energy. Events which had exactly two such cones were selected as e"'"e~ r"*"r~ 
candidates. The direction of each r was approximated by that of the total momentum 
vector of its cone of particles. Events were accepted if the average value of | cos^| for the 
two r jets, |cos6'av|, satisfied |cos6'av| < 0.85. 

Most of the remaining background from two-photon processes was rejected by a cut on 
the acollinearity and acoplanarity angles of the two r cones: the acoUinearity angle, in 
degrees, was required to satisfy 

^acoi < (180° - 2tan"^(2mzv^/(s - ml))) + 10° 

and the acoplanarity angle was required to be less than 30°. The value of the cut on 
acollinearity was chosen such as to include the peak from radiative return events at each 
energy; it is 50° at 130 GeV rising to 78° at 172 GeV. 

Remaining background from e"'"e~ — >■ e"'"e~(7) events was removed by rejecting events if 
the ratio of the electromagnetic energy to the track momentum in both of the r cones 
was between 0.9 and 1.1, as expected for an electron. 
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Finally, at 161 and 172 GeV, events classified as W-pair candidates according to the 
criteria in reference |^0| were rejected. 



The effective centre-of-mass energy of the e"'"e~ collision was estimated from the directions 
of the two r jets, as described for e"'"e~ events in section p73| . The distribution for the 172 GeV 
events is shown in figure §(d). A non-radiative sample of r^r~ events was selected from the 
inclusive sample by requiring s' / s > 0.8. 

The numbers of events selected at each energy, together with the efficiencies and feedthrough 
of events from lower s' into the s' / s > 0.8 samples, all determined from Monte Carlo simulations, 
are shown in table ^ 

The remaining background, which amounts to 5-13% in the inclusive samples and 2-7% in 
the non-radiative samples, is mainly from two-photon interactions; there are also contributions 
from electron and muon pairs. The total background contributions are shown in table |^, 
together with numbers of selected events and resulting cross-sections. 

The main sources of systematic uncertainty in the cross-section measurements arise from 
the efficiency and background estimation. The error in the efficiency has been estimated using 
high statistics samples of LEPl data, that in the background by comparing data and Monte 
Carlo distributions of the selection variables after loosening some of the cuts. 

The observed angular distribution of the r~ is shown in figure ^(c) for the s'/s > 0.01 
sample and figure |^(d) for the s'/s > 0.8 sample, for all centre-of-mass energies combined. 
Monte Carlo events were used to correct for efficiency and background, including feedthrough 
of events from lower s'/ s into the non-radiative samples. The corrected angular distributions 
at each energy are given in table |^. The forward-backward asymmetries were evaluated by 
counting the corrected numbers of events, as for the muons. Systematic errors were assessed 
by comparing different methods of determining the asymmetry: using tracks, electromagnetic 
clusters or both to determine the r angles. The total systematic error, including the contribution 
from the correction for interference between initial- and final-state radiation, is below 0.01 in 
all cases, much smaller than the statistical errors. The measured values are shown in table ||. 
In the same way as for the muons, we combine the 130 and 136 GeV data for the asymmetry 
measurements. From table |^ it can be seen that for the non-radiative sample at 133 GeV there 
is only one event in the backward hemisphere which after correction for efficiency, background 
and acceptance would yield an unphysical value of the asymmetry. 

Combined asymmetries from the fi^fi^ and r+r^ channels were obtained, assuming /x-r 
universality, by forming a weighted average of the corrected numbers of forward and backward 
events observed in the two channels at each energy. The combined values are shown in table 

3.6 The Fraction of bb Events 

To measure i?b, the ratio of the cross-section for bb production to the hadronic cross-section, we 
have performed b-fiavour tagging for the hadronic events with s'/s > 0.8, selected as described 
in section ^]2|. In addition we rquire at least seven tracks that pass standard track quality 
requirements, and the polar angle of the thrust direction to fulfil | cos6\ < 0.9 for the 161 and 
172 GeV data, |cos0| < 0.8 for the 130-136 GeV data. This acceptance cut ensures that a 
large proportion of tracks are within the acceptance of the silicon microvertex detector, which 
had a different geometry for the two sets of data. 
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Table 9: Numbers of selected events, efficiencies, backgrounds, feedthrough of events from 
lower s' to the s'/s > 0.8 samples and measured cross-sections for t~^t~ events. The errors on 
efficiencies and background include Monte Carlo statistics and systematic effects. The cross- 
sections labelled apneas are the measured values without the correction for interference between 
initial- and final-state radiation, those labelled cXcorr are with this correction. In some cases, 
the results are the same to the quoted precision. The ffist error on each measured cross- 
section is statistical, the second systematic. The Standard Model predictions, (Jsm, are from 
the ZFITTER ||25| program. 



The b-tagging technique is based on the relatively long lifetime (~ 1.5 ps) of bottom hadrons, 
which allows the detection of secondary vertices significantly separated from the primary ver- 
tex. The primary vertex for each event was reconstructed using a minimization method 
incorporating as a constraint the average beam spot position, determined from tracks and the 
LEP beam orbit measurement system. Although the beam spot is less precisely determined 
than at LEP 1, the resulting error on the primary vertex position is still small compared to the 
errors on the reconstructed secondary vertex positions. The secondary vertex reconstruction 
was the same as adopted in but the minimum number of tracks forming a vertex was 



reduced from four to three. Vertices were reconstructed in the x-y plane. Tracks used for 
secondary vertex reconstruction were required to have a momentum greater than 500 MeV. In 
addition, the impact parameter in the x-y plane relative to the reconstructed primary vertex 
was required to satisfy |(io| < 0.3 cm, and its error a^o < 0.1cm. This mainly removes badly 
measured tracks and, for example, tracks from K° or A decays. 

For each reconstructed secondary vertex, the decay length L was defined as the distance 
of the secondary vertex from the primary vertex in the plane transverse to the beam direc- 
tion, constrained by the direction of the total momentum vector of the tracks assigned to the 
secondary vertex. The decay length was taken to be positive if the secondary vertex was dis- 
placed from the primary vertex in the same hemisphere as the momentum sum of the charged 
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particles at the vertex, and negative otherwise. The distribution of decay length significance, 
defined as L divided by its error a^, combining data from all centre-of-mass energies, is shown 
in figure ^(a), superimposed on the Monte Carlo simulation. 

A 'folded tag' was used in this analysis in order to reduce the light fiavour component 
and the sensitivity to detector resolution uncertainties. Each hadronic event is divided into two 
hemispheres by the plane perpendicular to the thrust axis, and the hemispheres are examined 
separately. Each hemisphere is assigned a 'forward tag' if it contains a secondary vertex with 
a decay length significance L/ct/^ > 3, or a 'backward tag' if it contains a vertex with a decay 
length significance L/ctl < —3. Neglecting background in the hadronic sample, the difference 
between the number of forward and backward tagged hemispheres Nt — in a sample of Ahad 
hadronic events can be expressed as: 

Nt-Nt = 2N^^d [(eb - eb)i?b + (cc - e,)R, + (e^. - Cud.) (1 - Rh - Rc)] (5) 

where (eb — Cb), (cc ~ ^c) and (euds — ^uds) are the differences between the forward and backward 
tagging efficiencies. The difference for bb events (eb — eb) is about a factor of five bigger 
than that for cc events (ec — Cc), and a factor of fifty bigger than that for light quark events 
(cuds — Cuds) • -Rc is the ratio of the cross-section for cc production to the hadronic cross-section 
and was computed using ZFITTER. Due to the limited statistics compared with the LEPl 
data, a double tag technique cannot be applied in this analysis and one has to rely on a single 
tag method. Therefore the hemisphere tagging efficiency differences were determined from 



Monte Carlo, and are shown in table |I0|; the efficiencies vary only slightly with energy. The 
errors on these are predominantly systematic. The largest contributions to the systematic 
errors come from the uncertainties in Monte Carlo modelling of b and c fragmentation and 
decay, and from track parameter resolution. The b and c fragmentation and decay parameters 



were estimated by following the prescriptions of reference [42]. The effect of track parameter 



resolution was evaluated by varying the resolution in the transverse plane by 20%, in analogy 



to the procedure described in reference At the three centre-of-mass energy points the 

expected contribution from four-fermion background was subtracted, as described above for 
hadronic events. Within this background, only W-pair events are expected to contribute to the 
tagged sample. The probability for a W-pair event to be tagged was estimated from Monte Carlo 
to be (7.8±0.4)% at 161 GeV and (8.3±0.4)% at 172 GeV. The errors reflect the uncertainty 
of charm fragmentation in the W hadronic decay. After four-fermion background subtraction, 
b purities of the tagged sample of the order of 70% are obtained. 

The numbers of selected events, tagged hemispheres and resulting values of Rh are shown^ in 



table |TU|. The systematic error on i?b is dominated by the uncertainty on the tagging efficiencies. 
The other important systematic contributions result from Monte Carlo statistics and detector 
resolution. To check the understanding of the systematic errors, the analysis was repeated on 
data collected at the Z peak during 1996. The resulting measurement of R^, is in excellent 



agreement with the OPAL published value |4l|], differing by (0.7±1.7)%, where the error is 
purely statistical. 

The measured values of R^ at each energy are compared to the Standard Model prediction 
in flgure ^(b). Values for the bb cross-section, derived from the measurements of the hadronic 
cross-section and Rh, are given in table 

^The results presented here supersede those in reference [||. In particular the statistical errors on 
reference H were calculated incorrectly. 
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bb 




133.17 GeV 


161.34 GeV 


\12.V1 GeV 


Events 


255 


328 


296 


Forward tags 


61 


76 


66 


Backward tags 


10 


20 


25 


(eb - eb) 


0.414±0.023 


0.402±0.021 


0.395±0.022 




n n75+n nnfi 


W.W 1 i7_l_W.WWW 


n 080+0 006 

w . WO w_i_ w . www 


(euds e^ds) 


0.0059±0.0015 


0.0074ib0.0019 


0.0085±0.0021 


Dmeas 


0.199±0.040±0.013 


0.168±0.040±0.011 


0.136±0.048±0.010 


Dcorr 


0.195±0.039±0.013 


0.162±0.039±0.011 


0.131±0.046±0.010 


pSM 


0.182 


0.169 


0.165 


^hh (pb) 


12.5±2.6±0.9 


5.8±1.4±0.4 


3.5±1.4±0.3 


-bT (pb) 


12.7 


5.7 


4.6 



Table 10: Numbers of selected events, forward and backward tags, tagging efficiency differences 
and measured values of i?b- The values labelled B^^^ have not been corrected for interference 
between initial- and final-state radiation, those labelled R^" include this correction. The value 
of the bb cross-section, after correction for interference, is also given. The errors on the tagging 
efficiency differences include Monte Carlo statistics and systematic effects, the latter being 
dominant. The first error on or is statistical, the second systematic. The Standard 
Model predictions, -R^^, o"^^, are from the ZFITTER p5[ program. 



4 Comparison with Standard Model Predictions 



We compare our measurements with Standard Model predictions taken from the ALIBABA 
program for electron pairs, and the ZFITTER program for all other final states, with in- 
put parameters mz=91.1863 GeV, mtop=175 GeV, mHiggs=300 GeV, aem('^z)=l/128.896 and 
«s(^z)=0.118. We use ZFITTER version 5.0, with a small modification to the code to ensure 
that the s' cut is applied to fermion pair emission in the same way as it is to photon emis- 



sionQ. Measured values of cross-sections, presented in tables 0, ^, |^, |] and |T0|, are shown in 
figure |10[ The measurements are consistent with the Standard Model expectations. The asym- 
metry measurements, presented in tables ^ and ^, are shown in figure 11, while the corrected 



angular distributions for hadrons are shown in figure |I2| and for electron pairs in figure |I3. 
We have combined the differential cross-sections for muon and tau pairs, and show the average 
in figure O. The measured angular distributions and asymmetry values are in satisfactory 



agreement with the Standard Model expectations. 



In figure |I5| we show -Rinc, defined as the ratio of measured hadronic cross-section to the 
theoretical muon pair cross-section, as a function of centre-of-mass energy. The muon pair 
cross-sections are calculated using ZFITTER, as described above. The hadronic cross-sections 
used here are somewhat different from the measurements presented in section p]^. We use 



^We use default flag settings, except B0XD=1, C0NV=1, INTF=0 and FINR=0. The effect of the BDXD flag is 
to include the contribution of box dia gram s, which are significant at LEP2 energies; the setting of the other 
flags was discussed in sections Ifl and 2.3 . 
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an inclusive cross-section, (j(qqX), which is measured in a similar manner to the inclusive 
hadronic cross-section described above, but without subtraction of the W-pair contribution. 
The observed cross-section is corrected using an efficiency which includes the effect of W-pair 
events to give a total cross-section which is the sum of the two-fermion cross-section plus the 
cross-section for W-pair production with at least one of the W bosons decaying hadronically. 
This cross-section is thus an inclusive measurement of hadron production in e"'"e~ annihilation, 
in which production thresholds (e.g. for WW, ZZ or new particles) can be seen. The measured 
values of this cross-section and the ratio i?inc are given in table |ll]. In figure |15| the measured 
values of Rinc are compared with the prediction of ZFITTER, which does not include W-pair 
production, and also with a theoretical prediction Rf^ including the expected contributions 
from WW and ZZ events, calculated using GENTLE [45| and FERMISV [14| respectively. The 
effect of W-pair production is clear. 

In figure |T5| and table |TT] we also show -RBorn, the ratio of the measured hadronic cross-section 
for s'/s > 0.8, corrected to the Born level, to the theoretical muon pair cross-section at the 
Born level. The correction of the measured cross-section is performed using ZFITTER, and for 
both the numerator and denominator 'Born level' means the improved Born approximation of 
ZFITTER. Far below the Z resonance, this ratio becomes the usual that has been measured 
by many experiments at lower energy. Some of these low energy measurements are also 
shown in figure The measurements close to the Z peak have been corrected to our definition 
of -RBorn for this figurc. 



qqX 




130.25 GeV 


136.22 GeV 


161.34 GeV 


172.12 GeV 


O-corr (pb) 

Rinc 
pSM 


317±11±5 
14.4±0.5 
15.0 


264±10±4 
14.0±0.6 
14.5 


153±4±2 
13.6±0.4 
13.6 


138±4±2 
14.4±0.5 
14.2 


-RBorn 

pSM 

-'I'Born 


7.9±0.6 
9.5 


9.0±0.7 
8.9 


7.9±0.5 
7.5 


6.9±0.5 
7.2 



Table 11: Measured cross-sections for qqX, after correcting for interference between initial- and 
final-state radiation. The ffist error shown is statistical, the second systematic. Also listed are 
the ratios i?inc and -RBorn, as defined in the text, where statistical and systematic errors have 
been combined, and their respective Standard Model predictions. 



4.1 Influence on Electroweak Precision Measurements 

In references [|l],0 non-radiative data above the Z resonance were used to constrain the size 
of the interference terms between photon exchange and Z exchange processes, which have 
amplitudes of similar magnitude. Using the ZFITTER |^ and SMATASY programs, we 
have repeated the model- independent fits to OPAL data described in references [|],||, including 
all the measurements of the non-radiative hadronic cross-section and combined and r+r~ 

asymmetry presented here. In addition, we have included the measurements of the muon and 
tau pair non-radiative cross-sections; these were not used in the fits in references [0,0]. In the fit, 
the parameters ^ and j^, determining the sizes of the hadronic and leptonic 7Z-interference 
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respectively, have been left free (see reference for more discussion of these parameters and 
details of the fit). In the Standard Model, and jf' have the values 0.216 ± 0.011 and 
0.799 ± 0.001 respectively, for a top quark mass of 175 GeV and Higgs boson mass of 300 GeV, 
where the uncertainties in this prediction come from varying mniggs in the range 70-1000 GeV, 
as{mz) in the range 0.112-0.124, andmz 



^top; 

,2 



and tte 



(mz) in ranges taken from reference PS 

The results of the fit, which has a of 81-5 for 112 degrees of freedom, are given in 
table |l^. For comparison, the table also shows the results of the fit presented in [|[| to OPAL 
data collected at LEP 1 alone, and the Standard Model predictions for ^ and jf'. Since the 
7Z-interference vanishes on the Z peak, the inclusion of data far away from the Z resonance 
considerably reduces the uncertainty of [||]. The inclusion of the 130-161 GeV data reduced 
the uncertainty on by 55%; a further improvement of about 16% is observed by including 
the 172 GeV data presented here. As shown in reference this improvement is much larger 
than that which would be obtained by the inclusion of the full LEP 1 off-peak data. The high 
energy data also reduce the correlation between fitted values of j^" ^ and the Z mass, as can be 
seen in table O and figure ITB. 



OPAL data sample 


•tot 
J had 


mz 
(GeV) 


;tot ^ 
J had' 

correlation 


•fb 

k 


LEPl (1989-92) |6,43,4|] 
LEP 1 + these data 


-0.18 ±0.68 

-0.02 ±0.26 


91.187 ±0.013 
91.184 ±0.009 


-0.70 
-0.31 


0.684 ± 0.053 
0.715 ± 0.043 


Standard Model 


0.22±0.01 






0.799±0.001 



Table 12: Fitted values of the hadronic 7Z-interference parameter, the Z mass, mz, and 
the leptonic 7Z-interference parameter, j^, using different OPAL data samples. The mz values 
are quoted for the s-dependent Z-width. 



4.2 Energy Dependence of oiem 

The cross-section and asymmetry measurements for s'/ s > 0.8 have been used to investigate the 
energy dependence of the electromagnetic coupling within the framework of the Standard 
Model. Here the main information does not come from the 7-Z interference, studied in the fit 
described above, but from the pure s-channel photon exchange in the leptonic cross-sections. 

For final state muon or tau pairs with s' ~ s, the photon exchange process, which is pro- 
portional to «^j^(Q^), dominates the Z exchange by factors of 2-4 at centre-of-mass energies of 
130-172 GeV. The leptonic forward-backward asymmetries, reflecting the 7-Z interference, are 
somewhat less sensitive since they depend only linearly on aem{Q'^)- For hadronic final states 
at s' ~ s the ratio between Z and photon exchange is about a factor of five larger than that 
for leptons, so that they are dominated by Z exchange, even at the highest energies presented 
here. In contrast to the charged leptons, quarks have a sizeable vector coupling to the Z, which 
depends on sin^6'w(Q^), which in turn is closely related to aem{Q^) via the Fermi constant 
Gp. For hadrons, changes to affect the photon amplitude and Z amplitude (via the vector 
coupling) in opposite directions. At 133 GeV, the hadronic cross-section is still sensitive to 
via the Z amplitude, while at larger centre-of-mass energies its sensitivity is considerably 
reduced due to the relatively increased photon amplitude. 

To fit for aem{y/s) at each centre-of mass energy, we form the between our measurements 
and the Standard Model predictions calculated as a function of aem(\/s) using ZFITTER, 
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keeping all other ZFITTER input parameters fixed to the values given above. The dependence 
of the vector coupling on is then effectively included in ZFITTER. In addition to performing 
a fit for each centre-of-mass energy we also perform a combined fit to all energies, in which 
runs with energy with a slope which corresponds to its fitted value. In this case, we quote as 
the result the value of at a centre-of-mass energy ^/s= 157.42 GeV, corresponding to the 
luminosity weighted average of 1/s. 

As inputs to the fits, we use the measured values of hadronic, muon pair and tau pair cross- 
sections, i?b, and the combined muon and tau asymmetry values, for s' / s > 0.8, giving five 
measurements at each energy. The measurements at 130 and 136 GeV have been combined for 
this analysis. The correlation coefficients between all the measurements are given in table |13 . 
The most significant correlations are between measurements of the same quantity at different 
energies arising from common efficiency and background systematics. The error on luminosity 
measurements gives rise to correlations between cross-section measurements at the same energy, 
and also, through the systematic error, to small correlations between cross-sections at different 
energies. The errors on the correction for interference between initial- and final-state radiation 
are assumed to be fully correlated. 

The results of the fits, for the separate centre-of-mass energies and the combined fit, are 
shown in table |T^, with the resulting values of aem(v^) plotted in figure [1^. The values are 
consistent with expectation. The asymmetry in the errors arises because the dependence of the 
asymmetries on asm is linear whereas the cross-sections and -Rb have a quadratic dependence. In 
table |l^ we show, for each measurement, the normalized residual before and after the combined 
fit and an estimate of the sensitivity to the value of aem- The normalized residual is the 
difference between measured and predicted value divided by the measurement error, while the 
sensitivity estimate used here is the number of experimental standard deviations by which the 
prediction for the quantity changes if l/aem(157.42 GeV) is varied by ±1.0 from its Standard 
Model value. We see that the measurements of lepton cross-sections are the most sensitive 
to the value of aem, while the influence of the hadronic variables is considerably smaller, as 
explained above. From the combined flt, we obtain a value of 

l/aem(157.42 GeV) = 119.9l||, 

where the error arises from the errors on the measurements. The error due to uncertainties 
in the input parameters to ZFITTER is negligible in comparison, being at most 0.04 for a 
variation of mniggs in the range 70-1000 GeV. 

The fit described above uses measurements of cross-sections which depend on the measure- 
ment of luminosity, which itself assumes the Standard Model running of ttem from (Q^ = 0) 
to typically = (3.5 GeV)^, where l/ttem — 134. Therefore it measures running only from 
Q ~ 3.5 GeV onwards. To become independent of this assumption, we have repeated the 
combined fit replacing the cross- sect ions for hadrons, muon and tau pairs with the ratios 
cr(/i/i)/cr(qq) and a (rr) /a (qq). This is possible since, above the Z peak, leptons and hadrons 
have a very different sensitivity to aem, as explained above. The increased correlations between 
the input quantities have been taken into account. This fit yields a result of l/aem(157.42 GeV) 
= 119.9l5;4 ± 0.1, with a x^/d-o.f. of 7.5/11, very close to the result obtained above but with 
somewhat larger errors. The second error comes from varying ag(mz) in the range 0.112 to 
0.124. This measured value of 1/aem is 2.6 standard deviations below the low energy limit of 
137.0359979±0. 0000032 [Q, indicating the running of ttem from Q"^ = independent of the 
luminosity measurement. There has been only one previous measurement where the run- 
ning of at large has been measured dominantly from the photon propagator and not, as 
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133.17 GeV 


161.34 GeV 


172.12 GeV 




(T(qq) 


Rv> 




a{TT) 




(T(qq) 


Ry> 


fT(/i,//) 


a{TT) 


.4kb 


fT(qq) 


Ri> 




a{TT) 




o-(qq) 


1.000 


-0.004 


0.005 


0.003 


-0.003 


0.135 


-0.006 


-0.001 


0.000 


-0.003 


0.118 


-0.005 


-0.001 


-0.000 


-0.003 


i?b 




1.000 


0.001 


0.001 


0.001 


-0.004 


0.085 


0.001 


0.001 


0.001 


-0.003 


0.067 


0.001 


0.001 


0.002 








1.000 


0.003 


0.001 


-0.000 


0.001 


0.012 


0.001 


0.001 


-0.000 


0.001 


0.014 


0.001 


0.002 


a{TT) 








1.000 


0.001 


-0.000 


0.001 


0.001 


0.022 


0.001 


-0.000 


0.001 


0.001 


0.020 


0.001 












1.000 


-0.002 


0.002 


0.001 


0.001 


0.004 


-0.002 


0.001 


0.001 


0.001 


0.004 


o-(qq) 












1.000 


-0.006 


0.002 


0.002 


-0.003 


0.122 


-0.005 


-0.001 


0.000 


-0.003 


Rh 














1.000 


0.001 


0.001 


0.002 


-0.005 


0.057 


0.001 


0.001 


0.002 


















1.000 


0.003 


0.001 


-0.000 


0.001 


0.014 


0.001 


0.002 


C7(rr) 


















1.000 


0.001 


0.000 


0.001 


0.001 


0.025 


0.001 






















1.000 


-0.002 


0.001 


0.001 


0.001 


0.004 


cr(qq) 






















1.000 


-0.004 


0.001 


0.001 


-0.003 


























1.000 


0.001 


0.001 


0.002 




























1.000 


0.002 


0.002 


a{TT) 




























1.000 


0.001 
































1.000 



Table 13: Correlation coefficients between the various measured quantities (s'/s > 0.8) at each energy. 



on the Z peak |51|, from the Standard Model relation between the weak and electromagnetic 
coupling constants. Using a slope corresponding to the fitted value, the combined fit result cor- 
responds to a measurement of l/aem("^z) = 121. 4^4:9 ±0.1. This measurement is completely 
statistics limited, essentially independent of mmggs, and does not rely on assumptions about the 
running of at low Q^, which is the dominant uncertainty in the current value of l/aeui{tnz) 
of 128.90±0.09 |5^. If the error can be reduced by more than one order of magnitude with 
future statistics, a measurement such as this will lead to an improved knowledge of l/a;em(mz). 





Fit 


Standard Model 


(GeV) 




xVd.o.f. 


1/ Olcra 


xVd.o.f. 


133.17 
161.34 
172.12 


116.9+|-^ 


4.2/4 
3.4/4 
0.9/4 


128.3 
128.1 
128.0 


6.2/5 
5.3/5 
0.9/5 


157.42 


119.9^ 


10.6/14 


128.1 


12.9/15 



Table 14: Results of fits for acm- The first three rows show the fits to data at each centre-of- 
mass energy, the last row the combined fit to all energies. The Standard Model values of 1/aem, 
and the between the measurements and the Standard Model predictions are also given for 
comparison. 



V~s (GeV) 


Measurement 


X(SM) 


x(fit) 


Sensitivity 


133.17 


(T(qq) 


-1.60 


-1.29 


0.047 




Rh 


0.31 


0.42 


0.011 




f7(/i/i) 


1.85 


1.37 


0.045 




cr(rr) 


-0.35 


-0.79 


0.042 






0.13 


0.33 


0.021 


161.34 


CT(qq) 


0.84 


0.80 


0.006 






-0.16 


-0.02 


0.016 






0.09 


-0.59 


0.068 




(j{tt) 


1.74 


1.29 


0.046 






-1.24 


-1.03 


0.023 


172.12 


o-(qq) 


-0.33 


-0.44 


0.003 






-0.73 


-0.60 


0.014 






-0.39 


-1.08 


0.070 




cr(rr) 


0.07 


-0.46 


0.054 






-0.34 


-0.14 


0.022 



Table 15: Normalized residuals x and sensitivities for the agm fit to the combined data. For 
each measured quantity we give the residual before and after the fit, and the sensitivity, as 
defined in the text. 



5 Constraints on New Physics 
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5.1 Limits on Four-fermion Contact Interactions 

In this section we use our measurements of non-radiative cross-sections and angular distribu- 
tions at 130-172 GeV to place limits on possible four-fermion contact interactions. This analysis 
is similar to those presented previously Because the Standard Model cross-sections de- 

crease as 1/s, the sensitivity of the measurements to the contact interaction increases with 
centre-of-mass energy and so the inclusion of data at 172 GeV is expected to give improved 
limits. In addition, the analysis of the hadronic cross-section in references assumed that 
the contact interaction is flavour blind. Here the study is extended to the case where the new 
physics couples exclusively to one up-type quark or one down-type quark. 



5.1.1 Theoretical Expectation 

The Standard Model could be part of a more general theory characterized by an energy scale A. 
The consequences of the theory would be observed at energies well below A as a deviation from 
the Standard Model which could be described by an effective contact interaction, as depicted 



in figure [T^. In the context of composite models of leptons and quarks, the contact interaction 
is regarded as a remnant of the binding force between the substructure of fermions. If electrons 
were composite, such an effect would appear in Bhabha scattering (e"'"e~ e"'"e~). If the other 
leptons and quarks shared the same type of substructure, the contact interaction would exist 
also in the processes e"'"e~ —>■ fi'^ fj,~ ,e'^e~ —>■ t^t~ and e"'"e~ — *• qq. Alternatively, a four- 
fermion contact interaction could be a good description of deviations from the Standard Model 
due to the exchange of a new very heavy boson of mass mx if ^ \/s- More generally, 
the contact interaction is considered to be a convenient parametrization to describe possible 
deviations from the Standard Model which may be caused by some new physics. The concept 
of contact interactions with a universal energy scale A is also used in ep and pp collisions to 
search for substructure of quarks or new heavy particles coupling to quarks or gluons. 

In the framework of a contact interaction |^3[ the Standard Model Lagrangian for e"'"e~ — *• ff 
is extended by a term describing a new effective interaction with an unknown coupling constant 
g and an energy scale A: 

= ^^^^ E^%ie.7^e.][f,7,f,] (6) 

where 5 = 1 for e"'"e~ e^e~ and 5 = otherwise. Here eL(fL) and eR(fR) are chiral- 
ity projections of electron (fermion) spinors, and rjij describes the chiral structure of the 
interaction. The parameters 7]^ are free in these models, but typical values are between 
— 1 and +1, depending on the type of theory assumed H]. Here we consider the same 



set of models as in reference HI. In addition we consider three other models: LL + RR 



(^LL = V^K = ±1, Vlr = Vrl = 0) and LR + RL (?7ll = Vrr = 0, Vlr = Vrl = ±1) [0,11 
which are parity-conserving combinations consistent with the limits from atomic parity viola- 
tion experiments, and Odb {vll = ±1, '7rr = ±4, ?7lr = Vrl = ±2) |56|. 

Equation can be rewritten in terms of the parameter e = (g^/47r)/A^. In this case the 
differential cross-section can be expressed as 

^^ = asMis,t) + C',is,t)e + C^^{s,t)e^ . (7) 
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Here t = —s{l — cos6')/2 and 6 is the polar angle of the outgoing fermion with respect to the 
e~ beam direction. The term describes the interference between the Standard Model and 
the contact interaction, the term is the pure contact interaction contribution. Their exact 
form depends on the type of fermion in the final state and the particular model chosen, and is 
given, for example, in reference ||^[|. If the underlying process is the exchange of a new heavy 
scalar particle, X, then equation (D and equation are good approximations of the process 
so long as mx ^ ^/s. Limits on the energy scale of the new interaction are extracted assuming 
gV47r = 1. 

In order to compare the models with the data the lowest order expression given in equa- 
tion was modified to include electroweak and QCD effects, and experimental cuts and 
acceptances were taken into account. The Standard Model cross-section ctsm has been calcu- 
lated using ALIBABA for the e"'"e~ final state and ZFITTER for all other final states. Standard 
Model parameters were fixed to the values given in section |. The errors on these quantities 
are negligible compared to the statistical precision of the data. 

The dominant part of the electroweak corrections is due to initial-state radiation, which 
was taken into account by convolving the theoretical cross-section with the effects of photon 
radiation according to [0]. Initial-state radiation was calculated up to order in the leading 
log approximation with soft photon exponentiation, and the order a leading log final state QED 
correction was applied. Other electroweak corrections were taken into account by evaluating 
the cross-sections with the appropriate value of sin^ 9^. For the hadronic cross-section, QCD 
effects were taken into account by multiplying the electroweak-corrected cross-section by 5qcd = 
l + a,/7i + l.m{aj7iy. 



5.1.2 Analysis and Results 



We have fitted the data presented here on the angular distributions for the non-radiative 
e"'"e~ — >• e"'"e~, e"'"e~ i e"'"e~ — r+r^ processes, the non-radiative cross-sections for 

e"'"e~ — s> qq, and the measurements of R\,. In the case of the leptonic angular distributions, a 
maximum likelihood fit was used, where the total sample of candidate events was fitted with 
the theoretical prediction plus the background. The number of events predicted in each bin of 
cos 6 is given by 



Nr^'ie.r, cos( 



:i + r) 



ak{e)Ek{cos6) + a 



predjA; , 
bgd ' 



COS( 



(8) 



where ak{£) is the cross-section at the centre-of-mass energy point k, which is a function of the 
free parameter e, Ek{cos 9) is the correction factor for experimental efficiency, crbgd^''^(cos 9) is the 
predicted background cross-section, and is the integrated luminosity. The background cross- 
section is evaluated including efficiency and detector acceptances for the different background 
processes. The factor r is a scaling factor which takes uncertainties from systematic errors 
into account. It allows the overall normalisation to vary within bounds set by the appropriate 
systematic errors 0|. We derive 95% confidence limits from the values of e corresponding to 
a change in the likelihood of 1.92 with respect to the minimum. In the case of the hadronic 
and bb cross-sections, we used a fit to the measured values, incorporating the correlations 
between the measurements as for the ttem fit- In this case the 95% confidence limits correspond 



* Equation (2) in reference 
replaced by 2s. 



has a typographical error: the factor 4s on the left-hand side should be 
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to a change in of 3.84. In the case of the hadronic cross-section, we fit both for the case 
where the new physics couples equally to all five flavours and under the assumption that the 
new interaction couples only to one flavour. 



The results are shown in table |l^ and are illustrated graphically in figure |19|; the notation 
for the different models is identical to reference [Q. The values for bb are obtained by fitting 
the cross-sections for bb production, and there is no requirement on whether or not the new 
interaction couples to other flavours. By contrast, those for up-type quarks and down-type 
quarks are obtained by fitting the hadronic cross-sections assuming the new interaction couples 
only to one flavour. Most of the fitted values of e show no significant deviation from zero, 
indicating agreement with the Standard Model. However, the quadratic form of the cross- 
section (equation (^) can lead to two local minima, and in a few cases the local minimum near 
zero gives the higher or — log(likelihood). As an example, in figure [2^ we show the negative 



log likelihood curve for the fit to the e^e~ angular distributions, and the curve for the fit 
to the hadronic cross-sections assuming couplings to one up-type quark only, in both cases for 
the VV model. In the former case the positive interference between the Standard Model and 
contact interaction amplitudes results in a curve with only a single minimum, whereas in the 
latter case the negative interference results in two local minima. As can be seen in figure pO|(b), 
one of these is near zero, but the second one gives the lower x^. In this case, we quote 0.39 as 
the central value of e, while the limits on A are derived from the points e = 0.44 and e = —0.05 
where Ax^ is 3.84 above the right-hand and left-hand minimum respectively. 

The two sets of values A+ and A_ shown in table |1^ correspond to positive and negative 
values of e respectively, reflecting the two possible signs of r]ij in equation (j^). As before, the 
data are particularly sensitive to the VV and AA models; the combined data give limits on A 
in the range 4.7-7.7 TeV for these models, roughly 1 TeV higher than those for 130-161 GeV 
data alone. For the other models the limits generally lie in the range 2-5 TeV, approximately 
0.5 TeV above those from previous data. Those for the Odb model are larger (3-12 TeV) 
because the values of the t] parameters are larger for these models. The limits obtained here are 
complementary to those obtained in pp collisions ||5^. They are superior to published contact 



interaction bounds from ep collisions at HERA and severely constrain interpretations of 



the recently reported excess of events at high momentum transfer at HERA [Q. 



5.2 Limits on Heavy Particles 

The contact interaction analysis is an appropriate framework for searching for effects arising 
from the exchange of a new particle with mass mx 3> -^/i. For lower mass ranges, a/s ^"^x < A, 
we search for signs of new physics not within the framework of the contact interaction, but under 
the explicit assumption that the new phenomena are due to a heavy particle, which couples to 
leptons and quarks. Such a particle could be a leptoquark |Q or a squark in super symmetric 
theories with i?-parity violation [^. Beyond the kinematic limit for direct production, such a 
new particle might be seen through a change of the total cross-section in the process e"'"e~ qq 



via a t-channel exchange diagram as depicted in figure |T^. The allowed leptoquark states can be 
classified according to spin and weak isospin I. We denote scalar particles S/ and vector particles 
V/. Isomultiplets with different hypercharges are distinguished by a tilde, as in reference 0. 
The coupling strength of the leptoquark is given by a coupling constant gL or gR, where L and 
R refer to the chirality of the lepton. The two scalar states Sq and S1/2 and the two vector 
states Vo and V1/2 can in principle couple to both left and right chiralities. The product gL -gR 
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Table 16: Results of the contact interaction fits to the angular distributions for non-radiative e+e~ — > e''"e~, e"'"e~ — > /i"*"//", e+e~ — > 
T^r~ , the cross-sections for e"'"e~ — > qq and the measurements of i?b presented here. The combined results include all leptonic angular 
distributions and the hadronic cross-sections. The numbers in square brackets are the values of [^7ll!''7rr,''7lr, ''7rl] which define the 
models. Eq is the fitted value of £ = 1/A^, A± are the 95% confidence level limits. The units of A are TeV, those of Eq are TeV~^. 



is constrained very strongly by low energy processes 0. Therefore in this analysis only one 
coupling at a time is assumed to be non-zero. 

In a t-channel reaction the exchange of Sq and S1/2 are equivalent to the exchange of an 
i?-parity violating down-type squark and up-type antisquark respectively. The coupling to 
electrons is given by the term \[ji^LiQjD^ of the superpotential |6^, where the indices denote 
the family of the particles involved. Li and Qj are the SU(2) doublet lepton and quark super- 
fields and Df. denotes a down-type antiquark singlet superfield. The coupling constant A'^^^ is 
equivalent to gL in the leptoquark case since the lepton involved is left-handed. In the limit of 
very large masses, rrix 3> y/s, the leptoquark coupling constant gLQ is related to the contact 
interaction coupling constant gci by gLq/'^x = gci/^- 

The cross-section for e+e~ — qq including t-channel exchange of a leptoquark has been 
calculated by several authors, for example references p|-H. It can be written as 



J n 100 2-^ Pik\Jik\ 

dcost^ IzStts . ,j D 

where the colour factor A*",, = 3, pik are the spin density matrix elements and fik are the helicity 
amplitudes for the process e^e"*" — > q^q. The helicity amplitudes used here are taken from 
reference 0. 

In order to compare the model with data the lowest order expression given above was 



evaluated taking electroweak corrections into account as described in section |5Tl| for the contact 
interaction analysis. No QCD corrections to the heavy particle propagator or vertex corrections 
were included. Standard Model cross-sections were calculated using ZFITTER. 



5.2.1 Analysis and Results 

In order to derive a limit on the mass and coupling of a heavy particle which couples to quarks 
and leptons, we calculate the between the measured non-radiative hadronic cross-sections 
and the model predictions. The predicted cross-section depends on the mass of the heavy 
particle and its coupling constant. Varying the mass in steps of 25 GeV, for each value we 
find the coupling constant which minimizes and determine the 95% confidence limit on 
the coupling corresponding to a change in of 3.84. We assume the presence of only one 
leptoquark multiplet at a time. As mentioned above, only one coupling, either gL or gR, is 
assumed to be non-zero. In the case of Si, Vi with gL, and S1/2, V1/2 with g^, two states of 
the isospin multiplet contribute to the hadronic cross-section. Their masses are assumed to be 
degenerate. 

We perform the analysis for two cases. In the first case we consider all non-radiative hadronic 
events, assuming a non-vanishing coupling to one quark family. In the second case we use the 



cross-sections for bb production as described in section ^]6| considering all possible leptoquark 
couplings to the b quark. 



The results from the analysis of the hadronic cross-sections are presented in figure |21|(a) 
and (b) for the possible scalar leptoquark states. The 95% confidence hmits are shown as a 
function of the mass rrix and the coupling constant gL or g^ of the new particle. The limits on 
gL for So are equivalent to limits on the Yukawa couplings A'^^-;, of supersymmetric models with 
i?-parity violation with j = 1,2 and k = 1,2, 3, where k is the family index of the exchanged 
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squark. The limits on S1/2 derived from the hadronic event sample are also limits on X[jf^ with 
j = 1,2,3 and k = 1,2, where j is the family index of the exchanged squark. Results from the 
analysis of the hadronic cross-sections for the vector leptoquark states are shown in figure p2|(a) 
and (b). We do not show limits on the So (Vq) leptoquark with couphng gR (gL) because the 
effect of these particles on the hadronic cross-section at these energies is too small. 



Figures 21(c) and ^(c) give the limits on scalar and vector leptoquarks respectively, derived 
in the analysis of the bb cross-sections. Most of these limits are considerably more stringent 
than those from the total cross-section analysis. The limits are weaker for the V1/2 leptoquark 
with coupling gR and the Vi leptoquark with coupling gL because in these cases two leptoquarks 
of the isospin multiplet contribute to the hadronic cross-section, but only one of these couples 
to b quarks. The limit on S1/2 is equivalent to a limit on X[j^ with j = 1,2,3. 

As can be seen in the figures our analysis is sensitive to leptoquark masses much higher than 
the beam energy. Direct searches at the Tevatron can exclude scalar and vector leptoquarks 



with Yukawa couplings down to (9(10 ) up to masses of ~ 225 GeV |6^. Our limits extend 
this excluded region for large couplings. 

As mentioned above, in the limit of very large mass, the leptoquark and contact interaction 
couplings are related by ghci/mx = gci/A. Where the cross-section formula of a leptoquark 
exchange corresponds to one of the investigated contact interaction models we have checked 
the consistency of the two analyses by calculating limits on the leptoquark coupling at masses 
of several TeV. We found excellent agreement. 



5.3 Search for Charginos and Neutralinos Decaying into Light Gluinos 

In general, a comparison between the observed number of non-radiative hadronic events and 
the Standard Model prediction may be used to constrain possible new particle production 
whose signature is similar to that of the non-radiative hadronic events. The upper limit on the 
excess in the observed number of events over the Standard Model prediction was determined 
to be 49.9 events at 95% confidence level, summing over all centre-of-mass energies and taking 
into account the correlated systematic errors between different centre-of-mass energies. The 
numbers of observed non-radiative hadronic events are listed in table 0, and the Standard 
Model prediction was calculated using ZFITTER with input parameters given in section An 
example with a particular model of new particle production is presented here. 

In supersymmetric extensions of the Standard Model, the supersymmetric partners of the 
gluons are Majorana fermions called gluinos. Experiments at the TEVATRON have excluded 



a gluino of mass up to about 150 GeV, independent of squark mass if missing energy 

is a relevant signature. However, a relatively long-lived very light gluino with mass less than 
1 GeV might not be completely excluded by any experiment | |66| | . 



The supersymmetric partners of the weak gauge and Higgs bosons mix to form charginos 
and neutralinos. Chargino and neutralino pairs {xtxj ^ hj = 1? 2 and x'iXp hj = 1) 2, 3, 4) can 
be produced in e^e~ annihilation via s-channel virtual 7 or Z, or in t-channel slepton exchange. 
The cross-section for chargino pair production is typically a few picobarns at a/s = 172 GeV, 
although, if the sneutrino mass is less than 100 GeV, the cross-section can be significantly 
reduced by interference between the s-channel and t-channel diagrams. In the usual SUSY 
scenario, xf decays into a virtual W boson and a neutralino. The cross-section for neutralino 
pair production is typically a fraction of picobarn at ^/s = 172 GeV for each process. Normally 
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X° {i = 2,3,4) are assumed to decay into a virtual Z boson and a lighter neutralino. The 
lightest neutralino is usually expected to be the lightest supersymmetric particle and is stable, 
however this is unimportant for the present analysis. 

Searches for charginos and neutralinos at various experiments have been carried out with an 
implicit assumption that light gluinos do not exist, so that the decay of charginos or neutralinos 



leads to signatures of large missing energy |67-7C]. If light gluinos exist and squark masses 



are comparable to the W boson mass, charginos and neutralinos decay dominantly into qq'g 
(where the q and q' are the same flavour in the neutralino case) through a virtual squark by 
the strong interaction ||66|. For larger squark masses, rriq ~ 1.5mw, the squark propagator 



reduces the decay width into qq'g, and this becomes comparable to the W*x° or Z*x° moc 
Since the gluino hadronizes into many hadrons, the signature of such chargino or neutralino 
events is very similar to ordinary qq events or WW hadronic decay events. Therefore, as 
proposed in reference [^], we have used the hadronic cross-section measurements for s' /s > 0.8 
presented in section to place limits on any possible contribution from chargino or neutralino 



pair production in which gauginos decay into qq'g. This method has the advantage of being 
insensitive to details of gluino fragmentation which would cause a significant uncertainty in a 
search for an excess of multijet events. 

The SUSY partners of SU(3), SU(2) and U(l) gauge bosons of the Standard Model (gaugi- 
nos) are assumed to have the same mass at the Grand Unification mass scale. Masses of these 
gauginos at the weak scale are determined by the renormalization group equations. As a result, 
the ratios of the SU(3), SU(2) and U(l) gaugino masses at the weak scale (M3 : M2 : Mi) are 
proportional to the strength of the gauge couplings at the weak scale (03 : a2 ■ ai). Therefore, 
if the gluino, which is the SU(3) gaugino, is very light (< 1 GeV) at the weak scale, all the 
other gauginos are also very light. We assume M3 = M2 = Mi = 0. Then the chargino and 
neutralino pair cross-sections depend on three SUSY parameters: the Higgsino mixing param- 
eter fi, the ratio of the two vacuum expectation values tan/? and the sneutrino mass rrii). We 



have calculated the cross-sections using the formulae in reference |Tl|]. We assume that the 
sneutrino mass is larger than the current limit of 43 GeV ||72|| . 

To determine the limits, the total number of xtxj ih 3 = 1) 2) and Xi'Xj ih 3 = 2, 3, 4) events 
expected at each centre-of-mass energy was calculated from the cross-sections, the integrated 
luminosity and the detection efficiency. We assume that the branching fractions of xf,2 QQ'g 
and X234 QQS 100% and that the detection efficiency for chargino events with s'/s > 0.8 
is the same as that for qq events at each centre-of-mass energy]^. We make a similar assumption 
about the detection efficiency for neutralino pairs, except when m^o < 10 GeV, where the decay 
multiplicity could be low. We conservatively assume the efficiency is zero for X2X2 events if 
1^0 < 10 GeV. If m^o 

A2 X.2 



m^o < 10 GeV. If m^o < 3 GeV, the decay X2 ~^ 7Xi be significant, which could affect 



^The two body decays xf ~^ t'^^ or — > are in principle possible if i) is lighter than xf^ or Xa respectively. 
However, assuming q and i) masses are almost degenerate in the parameter space considered in this analysis, as 
is the case in the Minimal Supersymmetric Standard Model, xf ~* Qq' or Xi ~^ ^¥1 followed by q ^ qg will be 
the dominant decay. 

^"Note that no cuts have been applied to remove hadronic W-pair events. If such cuts had been made, this 
might have removed chargino pairs. Therefore the assumption that the efficiency is the same is justified. A small 
decrease of 0(5 — 10%) on the s' cut efficiency is e.g. predicted in Monte Carlo simulations of WW^ K^K^+X, 
for high-mass K^K^ pairs, when the energy deposits in the hadron calorimeter are ignored. Assuming similar 
fragmentation behaviour for gluinos and strange quarks p6t , the missing energy caused by possible invisible 
decay products of light gluino bound states is therefore not expected to be large enough to reduce the efficiency 
of the s' cut significantly. 
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our detection efficiency; therefore we have conservatively assumed the efficiency to be zero for 
X2Xj events (j = 3,4) if m^o < 3 GeV. 

We ffist consider the case where the difference between the observed and expected num- 
bers of hadronic events is due to the production of charginos {xtxjjhj = 1)2) only. This 
corresponds to a SUSY model with the looser condition M3 = M2 = with no restriction on 
Ml 0]. The 95% confidence level lower limit obtained on the /i parameter as a function of 
tan j3 for this case is shown in figure EH. 



We then consider M3 = M2 = Mi = |^6[ and calculate the total number of chargino and 
neutralino events with the assumptions above. We derive limits on the SUSY parameter fi as 
a function of tan f3. The entire /z-tan f3 region is excluded. Therefore in this particular model 
these excluded parameter regions imply the exclusion of a light gluino under the assumptions 
that the branching fractions of xf,2 ~^ QQ'g X2 3 4 ^ are 100% and that the detection 
efficiency for chargino events and neutralino events (with m^o > 10 GeV) is the same as that 
for qq events with s'/s > 0.8 at each centre-of-mass energy. This limit would continue to hold 
even if smaller branching ratios or a slightly reduced efficiency of the s' cut were to decrease 
the total number of expected events by 23%. 



6 Conclusions 

We have presented new measurements of the production of events with two-fermion hadronic 
and leptonic final states in e"^e~ collisions at a centre-of-mass energy of 172 GeV, and updated 
similar measurements at 130-161 GeV. Special attention has been paid to the treatment of the 
interference between initial- and final-state radiation, and the contribution from four-fermion 
production. The measured rates and distributions are all consistent with the Standard Model 
expectations. In a model-independent fit to the Z lineshape, the inclusion of data at 172 GeV 
provides an improved constraint on the size of the interference between Z and photon ampli- 
tudes. Within the framework of the Standard Model, the data have been used to measure 
the electromagnetic coupling constant, giving l/aem(157.42 GeV) = 119.911;}, compatible with 
expectation. 

We have used these data to place limits on possible deviations from the Standard Model 
represented by effective four-fermion contact interactions. Limits are obtained on the energy 
scale A generally in the range 2-7 TeV, assuming g^/47r = 1. We have searched for the effect 
of a new heavy particle which might be exchanged in the t-channel. Limits are obtained on the 
coupling constants gL, gR between typically 0.2 and 0.6 for masses below ~ 200 GeV. These 
limits can be interpreted both as limits on leptoquarks, or in some cases as limits on squarks in 
super symmetric theories with i?-parity violation. Compared to previous searches we are able to 
improve existing limits in particular in the high mass region mx > 250 GeV. Limits have also 
been placed on chargino and neutralino decays to light gluinos in supersymmetric extensions 
of the Standard Model. 
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Figure 1: (a) The distribution of the difference in radial coordinate between the two clusters 
in Bhabha scattering events used for the sihcon-tungsten luminosity measurement at 172 GeV. 
Distributions of the radial coordinates of clusters for (b) the 'narrow' side and (c) the 'wide' 
side calorimeter. Distributions are shown after all cuts except the acoUinearity cut in (a) and 
the inner and outer radial acceptance cuts, on that side, in (b) and (c). Points show the data, 
while the histograms show the Monte Carlo expectation. The arrows show the positions of the 
cuts which define the acceptance. 
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Figure 2: The distributions of i?vis, the ratio of the total energy deposited in the electro- 
magnetic calorimeter to the centre-of-mass energy, and -Rbai, the energy imbalance along the 
beam direction, for hadronic events at each energy. Distributions are shown after all event 
selection cuts except the cut on that quantity. The positions of the cuts are indicated by the 
arrows. Points show the data and histograms the Monte Carlo expectations, normalized to the 
integrated luminosity of the data, with light shading indicating background from two-photon 
events and cross-hatching indicating other backgrounds, mainly four-fermion (including WW) 
and T~^T~ events. 
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Figure 3: Distributions of effective centre-of-mass energy \fs' reconstructed for hadronic events 
at 172 GeV. The distributions are shown (a) for all events, (b) for events with one photon along 
the beam axis, (c) for events with one photon in the electromagnetic calorimeter and (d) for 
events with more than one photon. In each case the points show the data and the histogram 
the Monte Carlo prediction normalized to the luminosity of the data. The arrows indicate the 
position of the cut used to select 'non-radiative' events. 
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Figure 4: (a) The distribution of the ratio of total electromagnetic calorimeter energy to 
the centre-of-mass energy for e+e~ —>■ e+e~ events in acceptance region B, |cos6'e-| < 0.7 
and 6'acoi < 10°, at 172 GeV. (b) The same distribution for the large acceptance selection, C, 
I cos6'| < 0.96 and 6'acoi < 10°, at 172 GeV. Distributions are shown after all cuts except the one 
on total electromagnetic calorimeter energy; the arrows indicate the positions of the cuts on this 
quantity, (c) The acollinearity angle distribution for events satisfying the inclusive selection, in 
the acceptance region A, | cos 6*2+1 < 0.9, | cos6'e-| < 0.9, at 172 GeV. The cut on acollinearity 
angle has not been applied, (d) Observed distribution of cos 9 of the outgoing electron in e+e~ 
events with ^acoi < 10°. AH centre-of-mass energies have been summed for this plot. In each 
case, the points show the data and the histograms the Monte Carlo expectations, normalized 
to the integrated luminosity of the data, with the background contributions shaded. 
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Figure 5: The distributions of y/s' for (a) hadronic events, (b) electron pair events with 
I cos 6*6+ 1 < 0.9, I cos 6*6- I < 0.9 and 6'acoi < 170°, (c) muon pair and (d) tau pair events at 
172 GeV. In each case, the points show the data and the histogram the Monte Carlo prediction, 
normalized to the integrated luminosity of the data, with the contribution from events with 
true s' / s > 0.8 shaded in (a), (c) and (d). The arrows in (a), (c) and (d) show the position of 
the cut used to select 'non-radiative' events. 
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Figure 6: (a) The time difference between hits in back;-to-back; TOF counters for muon pair 
candidates; data from all energies have been included in this figure, (b) The distribution of the 
ratio of the visible energy to the centre-of-mass energy for muon pair events at 130-136 GeV. 
(c) and (d) show similar distributions at 161 and 172 GeV. In each case, the points show the 
data and the histogram the Monte Carlo expectation, normalized to the integrated luminosity 
of the data, with the background contribution shaded. The arrows show the positions of the 
cuts, which for the i?vis distributions are positioned 0.15 below the expected value for radiative 
return events where the photon is undetected. 
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Figure 7: Observed distributions of cos 6 of the outgoing negative lepton in muon and tau pair 
events, (a) and (b) show muon pair events with s'/s > 0.01 and s'/s > 0.8 respectively, (c) 
and (d) tau pair events with s'/s > 0.01 and s'/s > 0.8. The points show data for all energies 
combined, while the histograms show the Monte Carlo expectations, formed by summing the 
predictions at each energy normalized to the measured integrated luminosity values. The 
background contributions (including feedthrough from lower s' in the s'/s > 0.8 cases) are 
shaded. 
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Figure 8: (a) The distribution of the ratio of the visible energy to the centre-of-mass energy for 
tau pair events. Note that visible energy here is defined as the scalar sum of track momenta plus 
electromagnetic calorimeter energy, with no correction for double counting, (b) Distribution of 
cos 9 of the missing momentum vector, calculated using electromagnetic clusters, for tau pair 
candidates. Data from all energies have been included in this figure. In each case, the points 
show the data and the histogram the Monte Carlo expectation, normalized to the integrated 
luminosity of the data, with the background contribution shaded. The arrows show the positions 
of the cuts. 

53 



-| — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

, OPAL 

-f 130-172 GeV 
m u,d,s,c WW 
u,d,s,c Z/y 
□ b 




0.25 



0.2 



(b) 



OPAL 



0.15 7 ^225 

o;22 



0.1 
0.05 




0.21^ 



\ 

0.21 



I ' I ' I ' I 



I . I . I . I 



88 90 92 94 

J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 

90 100 110 120 130 140 150 160 170 180 

V s/GeV 



Figure 9: (a) Decay length significance distribution for all centre-of-mass energies combined. 
The points show the data, the histogram the Monte Carlo prediction, (b) i?b as a function of the 



centre-of-mass energy. The points show the measurements presented here, and the value [41 



obtained on the Z peak. The errors are statistical and systematic, summed in quadrature. The 
solid line is the ZFITTER prediction for s'/s > 0.8. 
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Figure 10: Measured total cross-sections {s'/s > 0.01) for different final states at lower 
energies p6|, ^|4^ , and this analysis. Cross-section measurements for hadrons, bb, muon and 
tau pairs for s'/s > 0.8 from this analysis are also shown. The cross-sections for and r+r~ 

production have been reduced by a factor of ten for clarity. The curves show the predictions 
of ZFITTER for hadronic (solid), bb (dot-dashed), and r+r~ (dashed) final states and 

that of ALIBABA for the e+e" final state (dotted). 
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Figure 11: (a) Measured forward-backward asymmetry for electron pairs with | cos^^c-| < 0.7 
and ^^acoi < 10°, as a function of y^. The curve shows the prediction of ALIBABA, (b) Measured 
asymmetries for all (s'/s > 0.01) and non-radiative {s'/s > 0.8) samples as functions of y/s 
for /i^/i" and t~^t~ events (combined). The curves show ZFITTER predictions for s'/s > 0.01 
(solid) and s'/s > 0.8 (dotted), as well as the Born-level expectation without QED radiative 
effects (dashed). The expectation for s'/s > 0.8 hes extremely close to the Born curve, such 
that it appears indistinguishable on this plot. 
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Figure 12: Angular distributions for hadronic events with s' / s > 0.8. The data points show 
the measurements corrected to the direction of the primary quark/ ant iquark, and corrected to 
no interference between initial- and final-state radiation. The curves show the predictions of 
the ZFITTER program with no interference between initial- and final-state radiation (sofid) 
and with interference (dashed). 
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Figure 13: Angular distributions for c+c events with ^acoi < 10°. The points show the data, 
the curves the predictions of the ALIBABA program. 
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Figure 14: Angular distributions for fi^fi~ and t~^t~ events with s'/s > 0.8. The data 
points show the measurements, which have been formed as a weighted average of the /j^/j^ and 
r^r^ distributions corrected to no interference between initial- and final-state radiation. The 
curves show the predictions of the ZFITTER program with no interference between initial- and 
final-state radiation (solid) and with interference (dashed). 
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Figure 15: Ratio of measured hadronic cross-sections to theoretical muon pair cross-sections as 
a function of centre-of-mass energy. Values are shown for the inclusive cross-section, cr(qqX) 
and for the Born level cross-section, as described in the text. The dotted and dashed curves 
show the predictions of ZFITTER for these cross-sections, while the solid curve also includes 
the contributions from W-pairs calculated using GENTLE and from Z-pairs calculated using 
FERMISV. Measurements at lower energies are from references P^,B3|,H1,EB|. 
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Figure 16: Central values and one standard deviation contours (39% probability content) in the 
Jhad ^s. niz plane resulting from model-independent fits to the OPAL data samples described 
in section |4.1j . The horizontal band shows the Standard Model expectation j = 0.216±0.011 
for a top quark mass of 175±5 GeV and a Higgs mass range of 70-1000 GeV. 
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Figure 17: Fitted values of 1/aem as a function of Q, which is ^/s for the OPAL fits. The 
open circles show the results of the fits at each centre-of-mass energy, the closed circle the 
result of the combined fit. The value obtained by the TOPAZ experiment is also shown 
for comparison. The solid line shows the Standard Model expectation, with the thickness 
representing the uncertainty, while the value of l/aem(0) is shown by the dashed line. 
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Figure 18: Feynman diagram representing a generic contact interaction (left) and the exchange 
of a leptoquark or a squark in a t-channel diagram in an e+e^ colhsion. 
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Figure 19: 95% confidence level limits on the energy scale A resulting from the contact inter- 
action fits. For each channel, the bars from top to bottom indicate the results for models LL 
to Cdb in the order given in the key. 
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Figure 20: (a) Negative log likelihood curve for the contact interaction fit to the e+e~ angular 
distribution for the VV model, (b) curve for the contact interaction fit to the hadronic 
cross-section, assuming coupling to only one up-type quark, for the same model. 
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Limits on the coupling for Scalar Leptoquarks 
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Figure 21: 95% confidence exclusion limits on gL or gR as a function of mx, for various possible 
scalar leptoquarks. (a) and (b) show limits on leptoquarks coupling to a single quark family, 
derived from the hadronic cross-sections, (c) shows limits on leptoquarks coupling to b quarks 
only, derived from the bb cross-sections. The excluded regions are above the curves in all cases. 
The letter in parenthesis after the different leptoquark types indicates the chirality of the lepton 
involved in the interaction. The limits on the So and S1/2 leptoquarks can be interpreted as 
those on i?-parity violating (Ir and ui, squarks respectively. 
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Limits on the coupling for Vector Leptoquarks 
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Figure 22: 95% confidence exclusion limits on gL or gR as a function of mx, for various possible 
vector leptoquarks. (a) and (b) show limits on leptoquarks coupling to a single quark family, 
derived from the hadronic cross-sections, (c) shows limits on leptoquarks coupling to b quarks 
only, derived from the bb cross-sections. The excluded regions are above the curves in all cases. 
The letter in parenthesis after the different leptoquark types indicates the chirality of the lepton 
involved in the interaction. 
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Excluded regions of the tanP - ji plane for light gauginos 
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Figure 23: The 95% confidence level lower limit of as a function of tan /3 calculated assuming 
M2 = ^3 = 0, with no restriction on Mi. The region above the solid curve is excluded. 
The sneutrino mass is assumed to be larger than the current limit of 43 GeV, and the decay 
branching fraction of — > qq'g is assumed to be 100%. The contours for chargino masses 
m^± = 45 GeV, 50 GeV, 60 GeV and 70 GeV are also shown with dotted hues. With the 
additional restriction Mi — 0, the existence of light gluinos is excluded everywhere in the 
//-tan /3 plane. 
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